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A GENERALIZED DIGITAL CONTOURING PROGRAM 


By Ruben L. Jones 
Langley Research Center 

SUMMARY 

A generalized digital contouring program is presented and discussed. The con- 
touring program was developed by combining desirable characteristics from several 
existing contouring programs and can be easily adapted to many different research 
requirements. The overlaid structure of the program permits desired modifications to 
be made with ease. 

The contouring program performs both the task of generating a depth matrix from 
either randomly or regularly spaced surface heights and the task of contouring the data. 
Each element of the depth matrix is computed as a weighted mean of heights predicted at 
an element by planes tangent to the surface at neighboring control points. Each contour 
line is determined by its intercepts with the sides of geometrical figures formed by con- 
necting the various elements of the depth matrix with straight lines. 

Both input and output of the digital program are described, as well as the critical 
program variables and tests. The program variables, subroutines, overlaid programs, 
and the listing of each are described and cross-referenced. A sample problem composed 
of 552 data points is furnished. The resulting computer output listing, as well as the 
contour chart, is shown. 


INTRODUCTION 

Contour charts are usually thought of as being two-dimensional pictorial representa- 
tions of topographic formations of land masses relative to a smooth mathematical refer- 
ence surface (datum surface). Contour charts can also be useful in portraying data which 
are obtained during the course of research in various scientific disciplines and which 
would ordinarily be tabulated. Any set of data, such as barometric pressure relative to 
sea level, whether computed or observed, which can be referenced to a two-dimensional 
coordinate system can in principle be represented graphically in the form of a contour 
chart. Further, scalar variations in three-dimensional space are representable by a 
series of contour charts, each of which represents the observations within a section 
through the field. Mathematical formulations of complex problems can also be studied 
with the aid of contour charts. In fact, the utility of contour charts is limited only by the 



user's ingenuity in adapting the data to a form suitable for graphic representation 
(contouring). 

Since the usual technique for drawing contour charts requires that each scalar 
magnitude be plotted within a rectangular grid before contour lines are drawn, a digital 
contouring program is desirable from the standpoints of efficiency and manpower 
requirements. Several existing contourii^ programs (for example, ref. 1) were care- 
fully studied for accuracy, efficiency, and so forth, and each was found to have desirable 
characteristics. A digital contouring program incorporating the best characteristics of 
each as well as certain improvements in data processing was considered desirable. The 
objective was to develop a simplified, generalized digital contouring program, which 
could be readily adapted to many different research needs. This paper describes that 
program and the analytical concepts upon which it is based, as well as its operation. 

The program is written for the Control Data 6000 series computer systems and relies 
on library fimctions and routines found in references 2 and 3. 

SYMBOLS 

A,B,C,A',B',C',D' coefficients of the equation of a plane 
F function 

H elevation (scalar variation) above reference plane or line 

i,j variables 

l, j,k X-, y-, and z-component, respectively, of a vector 

k,k* last variable in a series 

m, Z dimensions of depth matrix 

N magnitude of a vector normal to surface or curve 

N vector normal to surface or curve 

P control point 

Q vector connecting two control points 
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R ratio between two radial distances, 

r radial distance between two points 

U discrepancy in a variable 

W computed weight 

X,Y rectangular coordinate axes with arbitrary origin 

x,y coordinates of a scalar variation in X-Y coordinate system 

Z fimction representing elevation (height) of a plane at an arbitrary point 

in that plane 

z height of plane at arbitrary point 

6 incremental change or discrepancy 

Subscripts: 

i,j,n arbitrary variables 

k,k' last variable in a series 

I particular variable or reference 

x,y X- and y-components of a vector 

DESCRIPTION OF PROBLEM AND METHOD OF SOLUTION 

The science of graphic representation of scalar variables as a function of two inde- 
pendent variables has been the result of two advances in computer technology — first, 
the development of supporting equipment which can interpret computer instructions to a 
mechanical plotter and second, the ability to process large masses of data at one time. 

Digital contouring programs, as a rule, perform two basic functions. First, a grid 
of equally spaced surface heights (terrain, for example) is generated from randomly 
spaced surface elevations. Since each coordinate of the grid point is equally distant from 
its neighbors, a grid for surface heights can be represented by a matrix of elevations 


3 



(the "depth matrix”). Thus, the value at each element (mesh point) of the depth matrix 
represents an elevation above or below a given reference, and the location of the element 
in the array (matrix) corresponds to its location within the grid relative to the grid 
origin (the (1,1) element of the matrix). The second function of a contouring program is 
plotting the contours of the surface represented by the depth matrix. 

Depth Matrix 

A segment of a surface profile within a section through a three-dimensional sur- 
face is shown in figure 1. Within the figure the control points Pj and are near 



Figiire 1.- Segment of a surface profile passing through a grid point and 

two control points. 


neighbors to the ith grid point, and the planes Zj and which are tangent to the 

surface at Pj and Pj+n> intersect the plane section (plane of the page) as shown. 

As can be seen, the heights of the planes at i are z- and and each is an approx- 

imation of the height of the profile at the grid point. Further, an improved approx- 
imation of is obtained from the average of the z values, and when each Z- is 
assigned a weight Wj, which is a normalized monotonic function of the distance r^ of j 
from the grid point i, a better estimate to the height of the profile at i is given by the 
expression 

k’ 

Hi = 2 WjZj (1) 

j=l 
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where k' represents both the total number of control points entering into the computa- 
tion of Hj and the control point farthest from i. If 



when 

r = ^(x - Xi)2 + (y - Yi)2 (2) 

then, since an expression for a normalized weight must satisfy the condition 
k’ 

j=i 

an expression for Wj, which is a monotonic function of the radial distance of j from 
the ith grid point, can be written as 



As can be shown, Wj diminishes from a maximum of 1 when j = i to a minimum of 
zero when j = k' . 

From the foregoing calculations, it is apparent that the height of a surface 
above a reference plane can be estimated as a weighted series of heights, each of which 
is the height predicted at the grid point by a plane tangent to the surface at each control 
point. The accuracy of each estimated height is a function of the number k' of neigh- 
boring control points used in the computation and of the distribution of the control points 
about the grid point, as well as their proximity to the grid point. Further, conditions 
were placed on the spacing of the control points relative to each other, and each height 
Hj is independent of grid spacing and orientation. 

As stated, the elevation at each grid point is estimated from a weighted series of 
heights, each of which is the predicted height of a tangent plane above a reference plane 
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(datum plane). If Fj^ is a function in x, y, and z which represents the equation of a 
plane tangent to the surface at the ith control point, then a general equation for the tan- 
gent plane is 


Fi(Xi,yi,Zi) = a: + b!x. + c|yj + d! z. = 0 

where zj is the height of the tangent plane above the datum plane at any • Fur- 

ther, since the surface gradient is proportional to a unit vector Nj normal to the sur- 
face at i. 


Gradient = VFj = b! i + Cj j + D- k ~ 
where V is a mathematical operator having the form 


ax ay az 


Thus, 





are the components of a vector normal to the surface at the ith control point. 

The particular equation of the tangent plane is determined at the ith control point 
by setting 


Al = § = - Vi - Vl - ^i 


( 4 ) 


where 


N . = B. = -i 

x,i 1 d! 

C- 


( 5 ) 
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and at any Xj,y^ 




!j(xj,y,) . -Aj - 


BjXj 


Vi = 


( 6 ) 


Each gradient for the ith control point is computed from a least- squares minimiza- 
tion process. If Qj is a vector joining a neighboring jth control point (a total of k) 
with the ith control point (the reference control point), the unit normal which mini- 
mizes the weighted sum of the squares of the projections of each Qj onto Nj is found. 
Analytically, 


^ ^j^^i " " Minimum 

j=l 


( 7 ) 


where Wj is the weight for each control-point pair as defined by equation (3) when 
k' = k. From equation (7), the I, j, and k components of are determined. 

The limiting indices k and k' imply that the number of control points utilized 
for the gradient computation and the grid-point computation may differ. In practice, a 
value k = k’ = 10 is usually satisfactory for all the computations in the expressions 
discussed above. It is, however, required that the distribution of control points about 
the ith point be reasonably uniform. For instance, if all the points should lie along a 
straight line or in a narrow band, the solution will be indeterminate. Thus, the digital 
contouring program searches the input data in an attempt to select only those neighboring 
points which will permit a good solution to each To search the data, neighboring 

points are selected by incrementing the inner and outer radii of a ring by one -tenth of 
the distance between grid points (the difference between the two radii). Further, to avoid 
the possibility of the control points being on a line with i, the sine of the angle included 
by each pair of vectors is not allowed to be less than 0.17. The inner and outer 
radii are incremented until either k control points have been selected or the outer 
radius has reached a value equal to 10 times the grid separation. If the distribution of 
points is poor, the gradient at the ith control point is rejected, and processing is con- 
tinued until each of the remaining gradients have been either evaluated or rejected. 

From the foregoing discussion, it is important to note that each of the expressions 
is independent of all assumptions pertaining to the depth-matrix orientation, origin, and 
so forth. Thus, the origin or spacing of the equally spaced surface grid can be changed 
at will without affecting the computed gradients of the represented surface. 

At each control point the digital computer program evaluates first the gradient as 
described and then each coefficient in equation (5) . The tangential-plane coefficients at 
each control point are then utilized in equation (1) to evaluate the scalar variation at each 
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mesh (grid) point in the depth matrix. All computations are performed in the units of 
the X- and y-coordinates for the control points. 


Numerical Techniques 

A basic assumption of the contouring routines in the program is that each element 
in the matrix is both positive and nonzero. Since a particular element evaluated as 
described can be legitimately either negative or zero, a positive bias is added to each 
element after evaluation. The bias is determined by the program as an integral multiple 
of 10 000. The integer multiplier is determined as the least value which when multiplied 
by 10 000 and added to either the minimum anticipated scalar variation or the actual min- 
imum data point (a program option) will result in a positive nonzero number. Unfilled 
elements in the depth matrix (a result of insufficient control-point density) can then be 
detected and filled by a subsequent routine. 

The depth matrix is subjected to numerous tests (for example, smoothness) and is 
corrected where desirable. If the relative difference in elevation between two consecu- 
tive elements of the depth matrix is found to exceed a maximum slope assigned by the 
user, the elements are adjusted relative to each other. The net effect of this operation 
is to smooth the depth matrix. Further, if elements within the depth matrix are found to 
be unfilled due to inadequate control-point distribution, the element is filled by use of 
either interpolation or extrapolation. Finally, dimensions of the matrix cannot exceed 
60 by 60 or be less than 3 by 3. 

One test is performed to aid the user in evaluating the quality of the matrix fit to 
the actual input data. The actual position of each control point is used to obtain an esti- 
mate from the depth matrix of its predicted scalar variation by interpolating between 
elements of the matrix. When the predicted and observed values are compared, a slight 
discrepancy should result, since the matrix is at best an approximation to the true 

surface. Each can then be categorized on the basis of whether it falls within the 

limits 
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( 8 ) 


where 


Uj = U._j + 6U (9) 

Thus, the distribution of the discrepancies is determined and plotted as a frequency dis- 
tribution curve. The program allows for all values of Uj where 



and 


6U = 20 


Contour Plots 

The program contours one 3 by 3 submatrix of the m by Z depth matrix at a time. 
Basically, by assuming no interpolation between mesh points, each 3 by 3 submatrix can 
be subdivided into a series of adjacent triangles by connecting each set of three elements 
with a straight line as shown in figure 2. Contour intercepts are then computed along 
each side of each triangle beginning with triangle I and progressing clockwise around the 
squares. In the figure, the order in which each side of a triangle is considered is also 
shown. Intercepts are determined by assuming that each side is sufficiently small to 
permit the elevation difference to be linearly related to the side length. (See ref. 1.) 



Figure 2.- The geometry of a 5 by 5 suhmatrix when subdivided by connecting 
groups of three elements with straight lines. 
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The procedure is always the same regardless of the number of elements in the 
depth matrix or the order of interpolation required between grid points. Interpolation 
between grid points is a linear expansion of each 3 by 3 submatrix into a larger sub- 
matrix. This process is repeated within the expanded submatrix until each element of 
the original 3 by 3 matrix has been considered. At this point a new 3 by 3 submatrix is 
processed in a similar manner, the third column of the first matrix becoming the first 
column of the second matrix. 

The depth matrix is processed three rows at a time. If the number of elements in 
any one direction of the depth matrix is even, the resulting contour chart is terminated 
prematurely in that direction. Thus, it is conceivable that the depth matrix may not be 
plotted in its entirety. 

PROGRAM ORGANIZATION AND DESCRIPTION 

The digital contouring program is written in FORTRAN IV for the Control Data 
6000 series computer systems, is overlaid, and requires approximately 65 000 octal 
words of storage. The overlay structure is shown in figure 3. 



Figure 5.- Overlay structure. 
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Each overlaid program and the subroutines associated with it are listed in table I. 
For each program or subroutine, the required storage (in octal words) and the fimction 
are given. Further, the page numbers serve as an index to the listing and the more 
detailed discussion of each nonsystem program and subroutine. System subroutines are 
not discussed in detail. 

Figure 4 is a flow chart showing the interrelation of the various overlaid programs 
and subroutines. The arrowheads indicate the direction of flow, and the circled numbers 
(1, 2, 3, etc.), with the exceptions of 8 and 9, indicate the order in which each overlay 
and subroutine is executed. The connectors 8 and 9 refer to disk files which were gen- 
erated in one overlay as input to another routine. The letters a, b, and c indicate the 
order of flow from and to common program connectors. 

The flow chart is rather detailed in its portrayal of the major functions of each 
program overlay and of each major subroutine. The various program options are indi- 
cated, as well as some of the more critical decisions. Emphasis has been placed on the 
program decisions which would terminate the execution of a routine in a normal fashion. 
It is necessary, however, to point out that certain program diagnostics which will abort 
the program are not shown but will be discussed subsequently. Additional flow charts 
are considered to be redundant and are not presented. 
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Figure 4 .- Generalized digital contouring program. 
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Program CONTOUR 


CONTOUR, the program executive, has three main functions. First, it reads the 
program options and constants. Second, it supervises the execution of each of the three 
lower level overlays. Third, storage for each array and program constant, which is 
generally common to each overlay, is set aside. 

OVERLAy ( LI NR t 0,0) 1 

PROGRAM CUN TOUR ( I NP UT , OU TP U T , T APE5= I NPUT , T AP E 6 = 0JT P J T , T APE T , T AP E 8 , 2 

1TAPE9) 3 

INTEGER FKMT 4 

DIMENSION FKMT(8) 5 

DIMENSION K.0RE ( <,) .KURNI A) 6 

REAL MINwRU2,MAXwRD2,MINwRD3,MAX<IRU3 7 

CCMMON /8LK2/ KOR E , KORN , I F y P , SOLE 8 

COMMON /3LK.3/ KM, LI 9 

COMMON /SLK4/ uISI,TANG,M 10 

COMMON /BLKi/ M I N WK D2 , M A Xini R 02 , M 1 NR RD 3 , M AX.JR U3 11 

COMMON /dLKb/ SP, SPP ,L AST , 01 SPGKD, LINEUNT 12 

COMMON /0LK7/ N , NN , P XI , P y 1 , X MAX , y M AX ,XM A y , y M Ay , C I 13 

CUMMiUN/riLKlA/ UNIT.CNVTOIN 14 

COMMON /aLKlS/ 1 R JT P , 1 ? T OP T , I SR TOP T 15 

CCMMON /BlKOOX ISRT 16 

COMMON /dLKAO/ KM A X , KM A X 1 , I 0 MX 17 

COMMI.N /LIMITS/ M 1 N 1 , M AX I , I A OD I 18 

COMMON /LlNKOC/ LINK, RECALL 19 

RLCALL=6LRECaLL 20 

LINK=4LL1NK 21 

3 FORMATIRAIO) 22 

KLAI) (5,3) FRMT 23 

READ(5,2) KMAX.KMAXi , lUMX 24 

2 FORMAT (214,110) 25 

READ (j.PrRMD ij I iT , UN 1 r ,C NIV I 01 1 , S PP ,C I , SC L E , U I SP GR 0 26 

l.TANG 27 

READ (5,4) ISR T, MINI , MAXI 28 

4 FORMAT ( 14 ,2 1 1 J ) 29 

L I NbuNi =01 SPGRU 30 

IPyp=liLASI=2 31 

READ (5,1) IkOT P , IP TuPT , I S3 T OPT 32 

1 FURMAT(3I4) 33 

CALL CVERLAy (L I NK , 1 ,0 , RECALL ) 34 

IF (ISRT.tJ.l) CALL CVERCAyl LINK.2 ,0 ,RECALL) 35 

CALL OVERLAy(L INK, 3, C , RECALL) 36 

STOP 37 

ENU 38 
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Program READIT 


READIT is the first overlay which is subordinate to CONTOUR. READIT is 
responsible for reading the X- and Y-position and scalar variation (Z-component) for 
each control point. As pointed out, the data may be read from tape or cards and may be 
either fixed or floating point. READIT applies the appropriate conversion to the input 
data. Further, the NW, NE, SW, and SE coordinates of the rectangle which will include 
all the control points, as well as IDIST (the distance of each control point from an arbi- 
trary origin), are determined by READIT. Each control point IX, lY, and IZ is then 
written on the appropriate file for use later along with IDIST and a code ICODE which 
serves to warn subsequent routines when the data have almost expired. 

OVlfKLAY (LiiNK.l.J) 39 

PKJGh<A'l KtAUir 40 

I.-JTbGtR FkMT 41 

Ul.-1c.Nii.lJiM f-R.'inS) 42 

JlMtiMSIuN KUkL( 4) ,K0RM( 4) ,/j;<D(3,4) 43 

rlF AL M A X W kU 2 t .'1 1 Jrt RU2 • M A Xw R 0 3 1 Nh KD 3 44 

CC.Yi'UJ.nI /ilLlM?/ KJKL ..xOkINi, IFYP .SOLE 45 

C JWMuK /oLK.4/ ul A ( 46 

Cn.MHUiM /ULin 5/ .-UM4NKJ2 ,.»lAXrt 3 J2,M1 ,MviKU3,MAXW.-<U3 47 

/aLKlSX IRUTP, IPrUPI , iSRTOPT 48 

CG-M.'-'ill \| /3LK30/ I SkT 49 

JUi-'.NuN YLI/ilTS/ .-iliNl .MAXI f lAJJl 50 

CUP.MJK /LIImKuJ/ LiMK, .kECAlL 51 

1 ruAMAf (8A10) 52 

RErtiNU 7 53 

ALPrlA= (-1 1 . 8/ laO .0 ) *( 3 . lAl 5926 54) 54 

ic;mt = 3 55 

IC0DE=0 56 

IF ( I kUTP. NL . 1 ) RLAU(5,1) FRMf 57 

100 CJ.MTIiMJE 58 

IE (ikJTP.Ew.l) oU TC lOCv. 59 

IE ( IPTJPI .Ob . 1 ) REAUIO.FEMD IX.IY.IZ 60 

IF ( IPTOPI.LF.l ) RLAi; (5 , FSMT ) X.Y.Z 61 

IE (cUF.a) 1002.1001 62 

1C 00 [ E I 1 prop r . Cb. 1 ) call <EC1.M(9.1.KK.IX.IY.IZ) 63 

IF t IPIGPl .LT .1 ) CALL lEClNi ( 9 . I . K. K . X , Y . Z ) 64 

lF(t(jF.9)lO0£i.lL01 65 

1001 CUIMTI.MJE 66 

IF I W IiM 1 .bw .a . A.MU .MAXI .£0 .0) 104.105 67 

1C5 CJ.mIINUE 68 

IF I Z.lT.i-UMI .'JR. Z.Gl.MAXI) 100. 104 69 

104 CJ.MTI.Nje 70 

IF ( 1 PT I . L I . 1 ) oJ TU 10u3 71 

X=FLL.ATiIX) 72 

Y^FLOAr ( I Y ) 73 

Z=FLU AT I i Z ) 74 

lOOOCUiMTIMJE 75 

IC,'MT= ICMI +1 76 

IF (IC.MT.bO.l) 300, 4bu 77 

3C0 CJimTINJE 78 

MINlAHU2=MAXrtR02 = Y 79 

HIl\rtKl)3 = MAX'AR03=X 80 


14 



400 CONTINUE 81 

ih (M IN WRD2 ,01. y) MINrti^J2 = y 82 

IF (MINWKD3.GT.X) MINwR03=X 83 

IF (MAXWRD2*LT. Y) MAXWRD2=Y 84 

IF ( MAXrtKl)3 .LT, X) MAXWRU3 = X 85 

IX^IFlX(X)$iy=IFlX(y)$IZ=IFlX(Z) 86 

X=X-l.OE20 87 

Y=Y-1.0E20 88 

IF ( ISKTOPT .LT. L) 1UIST=IFIX(Y/DIST) 89 

IF ( ISRTOPT.GE.l) 1 U IS T = I F I X ( X/D I S T ) 90 

CALL RECDUT ( 8 , 1 , 0 , 1 CODE , I X , I Y t 1 Z , ID 1ST) 91 

IF (ISkT.lt. 1) CALL KbCOUT ( 7 . I , 0 , I CODE , I X , I Y , I Z , I D I S T ) 92 

GO TO 100 93 

1002 CONTINUE 94 

KUREm=!<0RE(4) =MINWRD3 95 

KURE (2) =XGKE( 3) =MAXk^R03 96 

KORN( 1 )=KUKN(2 )=MAXWRD2 97 

KURN( 3) =K0RN(4) =MINWRD2 98 

IC0Dfc=9 99 

IUIST=0 100 

IX=IY=IZ=0 101 

IF (ISRT.LT.l) CALL KESGUT ( 7 , 1 , 0 f I C ODE , I X r I Y » I Z . I 0 1 ST ) 102 

CALL KECCUI ( 8 . i, 0, ICGOE, 1 X, I Y» IZf ID 1ST ) 103 

ICUDL=99 104 

IF (ISRT.LT.l) CALL KtCOUT ( 7 , 1 , 0 1 I C ODb , I X , I Y , I Z , I 0 1 S T ) 105 

CALL KECOUI ( 8 , 1 . 0 . I CJDb , i X , lY , IZ r I D IS T ) 106 

REWIND 8 107 

REWIND 7 108 

RETURN 109 

END 110 
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Program SORTIT 


SORTIT is an optional overlay which is responsible for setting up the appropriate 
input to SORT2, which in turn calls the SORT/MERGE routines of the computer system. 
These routines are very flexible and will utilize all the storage allotted to them. 

Regardless of which sorting option is chosen (fig. 4), ICODE is the primary vari- 
able. In this manner, the codes 9 and 99 are retained at the end of the data file. This 
fact in no way changes the indicated flow in the flow chart. 


OVERLAY lLIKKt2.0J 111 

PRGoRAH SORTIT 112 

COMMON /aCRla/ I KU T P , I P TOP T , I SKTOP T 113 

COMMON /LIKKOC/ LINK. RECALL 114 

UlPEiNSlON I SM(‘3» , 1FN(2) ,KEY( 16) 115 

REWIND 7 116 

REWiNU 8 117 

ISM(i)=l 118 

1SM(2)=A 119 

1SM(3)=60 120 

lSM(4)=lHf 121 

ISM(5»=1HB 122 

lE.NI 1 » = 5LTAPE7 123 

IFN(2)=6LTAPtd 124 

KEVm = lHA 125 

KEY(2)=1MX 126 

KEY(i)=2 127 

KEY(6)=1HX 128 

KEY(1Q)=1HX 129 

KEY(13)=1HA 130 

KEY(14)=1HX 131 

KEY(15)=5 132 

IF ( ISRTOPT.Ot. 1) GO Tu 1300 133 

KEY(5)=lh0 134 

KEY(7)=6 135 

IF I ISRTUPT.EO.-i) KEY(7)=4 136 

KE Y( 9) =1HA 137 

KEY(IiJ=3 138 

GO TO ICOl 139 

1000 CONTINUE 140 

KEY(5I=1HA 141 

KEY(7)=6 142 

IF (IS RTUPT.E0.2) KEY(7)=3 143 

KEY(9)=1HU 144 

KEY(11)=4 145 

1001 CONTINUE 146 

CALL S0RT2 ( ISM, IFN.KEY) 147 

END FILE. 7 148 

RETURN 149 

END 150 
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Program COUPLE 


COUPLE is the last of the three major overlays subordinate to CONTOUR. 

COUPLE is assigned the function of couplii^ the three lower level overlays. Further, 
storage for the depth matrix is set aside. 

OVERLAY (LINK, 3,0) 151 

PROGRAM COUPLE 152 

DIMENSION MAr(60,60) 153 

COMMON /BLKl/ MAT 154 

COMMON /LINKUC/ LINK, RECALL 155 

DO 1 1=1,60 156 

DO 1 J=l,60 157 

1 MArU,J)=0 158 

CALL OVERLAY(LINK, 3,1, RECALL) 159 

CALL OVERLAY ( L INK , 3 ,3 , REC ALL ) 160 

CALL OVERLAY(LINK,3,2,R£CALL) 161 

RETURN 162 

END 163 


Program MATCAL 

MATCAL furnishes a link between the subroutines GRIDIT and REGRID. 

OVERLAY (LINK, 3,1) 164 

PROGRAM MATCAL 165 

CCMMUN /LINKDC/ LINK, RECALL 166 

REWIND 7 167 

CALL GRIDIT 168 

CALL RbGRIO 169 

RETURN 170 

END 171 
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Subroutine GRIDIT 


GRIDIT is the first of several subroutines which have the responsibility for gen- 
eratii^ an adequate depth matrix. The indicated flow between connectors 5a and 5b 
(fig. 4) is descriptive of its functions. 

This routine is the largest consumer of computer time. The block in which the 
gradients are computed first searches the data for control points which are adequately 
distributed about the control point for which the gradient is desired. In the process of 
searching the data, no control point which has a relative difference in elevation suffi- 
ciently large to yield a slope greater than TANG (an input variable) is accepted. Further, 
IDMX is the maximum scalar difference allowed for the product of the distance R 
between the control points and TANG. Thus, it is possible for all the control points to be 
considered several times. To reduce this undesirable consumption of time, BIST can be 
increased or KMAX can be decreased. If the time problem is not solved by these actions, 
the data should be examined carefully. 

Each evaluated element of the depth matrix is made positive by adding a multiple 
of 10 000 to it. This constant is denoted as lADDI and is computed from the minimum 
scalar variation MINI. 

SUliRuUTINt GRIDIT 172 

DiyENSION RESULT (2J 173 

DIMENSION IKIIOO) 174 

DIMENSION RS(IOO) 175 

DIMENSION AA(3,3I ,88(3,1) , A( 1000 J ,B( 1000) ,C( 10001 ,0( 3, 1000) ,6SAVE( 176 

13, 100) ,USAVtI3) , IPIV0T( 3) , INDEX! 3,2) 177 

DIMENSION E(3) 178 

DIMENSION MAT(60,60) ,KUKN(4) ,KJkE(4) 179 

DIMENSION GRA0(3) ,GRADCJ(3,3) 180 

INTEGER S 181 

real MINwRU2,MAXWKD2,MlNWRD3,MAXrtR03 182 

COMMON /BLKl/ MAT 183 

COMMON /3LK2/ KURE , KUkN , I FTP , SCLE 184 

COMMON /8LR3/ JJ,KK 185 

COMMON /8LK4/ DIST.TANG 186 

COMMON /3LK5/ M I NWRD2 , M AXRRD2 , M INW RD3, M AXWR03 187 

COMMON /8LK6/ S P , SPP ,L AST , DI SPGRO 188 

COMMON /8LK40/ KMAX , KM AX 1 , 1 DMX 189 

COMMON /LIMITS/ M 1 N I , MAX I , I ADD I 190 

2000 FORMAT (7X#CORNER COORDINATES OF AREA TO BE CONTOURED*/) 191 

2002 FORMAT (* NORTH WEST CORNER — X=*E16.8* E16.8 ) 192 

2003 FORMAT I* SOUTH EAST CORNER — X=*E16.8* Y=* E16.8) 193 

2004 FORMAT (Ix, *Y DIMENSION OF DEPTH MATRIX IS* ,F 1 5. 5 , * INCHES* ) 194 

2005 FORMAT (IX, *X DIMENSION OF OEPTfl MATRIX IS* , F 1 5. 5 , * INCHES*/ ) 195 

2006 FORMAT ( IX , 14, *EL6MENTS IN Y DIRECTION OF DEPIH MATRIX*) 196 

2007 FORMAT (IX , 1 4 , *EuEMENTS IN X DIRECTION OF DEPTH MATRIX*/) 197 

ICNT=0 198 

RSAVE=DIST 199 

IDEPSUM=0 200 

PRINT 243 201 

240 FORMAT!* DEPTH MATRIX SIZE*/) 202 
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241 FURMAK/) 203 

242 FORMATi* THUS.4I 204 

LABL=0 205 

CSLOPE1=1.0/S3RT( l.0+TANG*#2» 206 

CSLOPE^O.9 207 

IF (CSL0PE.GT.CSL3PE1 ) CSL0PE=CSL0PE1 208 

243 FORMATllHlf* THE FOLLOWING IS INFORMATION PERTAINING TO THE DEPTH 209 

IMATRIX — S IZE. DIMENSIONS, ETC. *///) 210 

MINMAX=0 211 

IF (MINI. NE.O. OR. MAXI. NE.O) MINMAX=1 212 

IF (MINMAX.EO.O) 200,201 213 

200 MINI=200000 214 

MAXI=-20C000 215 

201 CONTINUE 216 

SP=OIST 217 

IPTOK’^O 218 

CALL CCNVERT(SP) 219 

PRINT 2000 220 

PRINT 2002 ,MINWR03 ,MAXWRD2 221 

PRINT 2003,MAXWRD3,M1NWR02 222 

CALL CONVERT (DiSPGRDI 223 

8 FORMAT (F5.2,F10.4,2I5, 15) 224 

10 FORMAT (2110) 225 

IF (IFYP) 4,3,4 226 

3 IY=KCRN(3»-K0RN(1 ) 227 

GO TO 2 228 

4 IY=K0RN(1)-K0RN(3) 229 

2 YMAX=FL0AT( IY» 230 

CALL CONVERT (YMAX) 231 

IX=KORE ( 3)-K0Rh ( 1 ) 232 

XMAX=FLOAT( IX) 233 

CALL CONVERT(XMAX) 234 

JJ=IFIX(YMAX/SP)+I 235 

KK=IF IX( XMAX/SP) +1 236 

IF (XMAX-FLUAT(KK-1)*SP.GT.0.0» KK=KK«-1 237 

IF (YMAX-FLOAT(JJ-1)*SP.GT.O.O) JJ=JJ+1 238 

IF (RK.LT.3.0K. JJ.LT.3) 301,302 239 

302 CONTINUE 240 

IF (KK.GI.oO.OR.JJ.GT .60) 230,231 241 

301 PRINT 300 242 

GU TO 304 243 

230 PRINT 232 244 

304 CONTINUE 245 

PRINT 2008 246 

PRINT 2004, YMAX 247 

PRINT 2005, XMAX 243 

PRINT 240 249 

PRINT 2006, JJ 250 

PRINT 2007, K'K 251 

232 FORMAH* SINCE UlSl IS TOO SMALL THE MAXIMUM DIMENSIONS OF THE DE 252 

IPTH MATRIX ARE EXCEEDED — 1(60,60)*) 253 

300 FORMAT)* SINCE DIST IS TOO LARGE THE MINIMUM DIMENSIONS OF THE DE 254 

IPTH MATRIX VIOLATED — 1(3,3)*) 255 

GO TO 215 256 

231 CONTINUE 257 

2008 FORMAT!//* DIMENSIONS OF REJUIRED PLOTTING SURFACE EXCLUSIVE OF SR 258 

IID*/) 259 

PRINT 2008 260 

DO 180 J=1,JJ 261 
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DO 180 K=1,KK 262 

MAT(J.K)=0 263 

180 CONTINUE 264 

K=1 265 

PRINT 2004, YMAX 266 

PRINT 2005, XMAX 267 

PRINT 2006, JJ 268 

PRINT 2007, KK 269 

TCST^l.OE+20 270 

S=-2 271 

DO 740 1=1,1000 272 

DO 740 J=l,3 273 

740 D(J,I)=TEST 274 

104 CALL RECIM 7,1 ,KKK, ICODE, IX , lY, IZ, IDISr ) 275 

IF (IFYP) 40,30,40 276 

30 IY=(-1)*IY 277 

40 CONTINUE 278 

IF UCODE.E0.99I GO TO 221 279 

IF I ICNT.EO.lOuO. AND. ICODE. lt. 9) 202,203 280 

202 PRINT 211 281 

GO lU 215 282 

211 FORMAT!* YOU HAVE OVER ICOO CONTROL POINTS.*/* ADDITIONAL CONTROL 283 

IPOINTS CAN NOT BE ACCEPTED.*/* STORED DATA MAY NOT ADEQUATELY REPR 284 

2ESENT AREA TO BE CUNTOURED.*) 285 

203 CONTINUE 286 

IF (LABL.FO.C) 204,205 287 

204 IF (MlNMAX.Eu.O) 212,213 288 

212 PRINT 206 289 

LA8L=1 290 

PRINT 242 291 

206 FORMAT!* YOU HAYE ELECTED TO REJECT NO CONTROL POINTS*) 292 

2C7 FORMAT!* YOU HAVE ELECTED TO REJECT ALL CONTROL POINTS WITH SCALAR 293 

1 VARIATIONS LESS THAN*I6* AND GREATER THAN*I6) 294 

208 FORMAT!* YOU HAVE A MINIMUM SCALAR VARIATION QF*I10) 295 

209 FORMAT!* YOU HAVE A MAXIMUM SCALAR VARIATION OF*llO) 296 

GO TO 21C 297 

213 PRINT 2C7 .MINI , MAXI 298 

LA8L=1 299 

PRINT 242 300 

IM1N1=20COOC$IMAXI=-2003DO 301 

210 CONTINUE 302 

214 FORMAT!* MUST TERMINATE PROGRAM EXECUTION WITH A MODE 4 FATAL ERRO 303 

1R*I 304 

GO TO 216 305 

215 CONTINUE 306 

PRINT 214 307 

AZER0=0.0 308 

AX=£!X/AZER!3 309 

AY=AX*AX 310 

216 CONT INUE 311 

205 CONTINUE 312 

IF IIC0DE.E0.9) 219,220 313 

219 IF IMINMAX.NE.O) PRINT 208,IMINI 314 

IF IMINMAX.NE.O) PRINT 209, IMAXI 315 

IF IMINMAX.EQ.O) PRINT 208, MINI 316 

IF IMINMAX.EQ.O) PRIM 209, MAXI 317 

IADDI = !-l)*!10uGO)*!MlNl/10000H-10000 318 

MIN1=M1N1 + I ADUl $MAXi=MAXH-I ADDI 319 

GO TO 221 320 
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220 IF ( MiNMAX.tO.O) 217,218 321 

217 IF IIZ.LT.MINI) MINI=IZ 322 

IF IIZ.GT.MAXI) MAXI=IZ 323 

GO TO 221 324 

218 CONTINUE 325 

IF (IZ.LT.IMINU IMIN1=IZ 326 

IF (IZ.GI.IMAXn IMAXI = IZ 327 

221 CONTINUE 328 

Bm=FLOAT(IX) 329 

B(2»=FL0AT( lY) 330 

IJUMP=-1 331 

IF (lCODE-9) 100,1133,105 332 

100 CONTINUE 333 

1133 CALL CKPUINT I 8 ( 1 ) , B I 2 ) , I Z , I PTUK , I CODE } 334 

IF ( IPTOK.EO.-l ) GO TO 104 335 

103 CONTINUE 336 

ICNT=1CNT+1 337 

1103 0 (l,ICNT) = 8m 338 

DI2, ICNT ) = 8I2) 339 

81 3»=FL0AT< IZ) 340 

D(3,1CNT)=B(3I 341 

IF (ICOOE.lt. 9) GO TO 104 342 

105 CONTINUE 343 

DO 700 J=l, ICNT 344 

IF ID(1, J) .EO.TESTI GO TO 700 345 

DO 750 11=1,3 346 

GRAO(II)=0.0 347 

88(11, 1)=0.0 348 

DO 750 I J=1 ,3 349 

AAdl.I J)=0.0 350 

750 CONTINUE 351 

DXMX=DYMX=-1000.0 352 

DXMN=DYMN=+1000.0 353 

DU 705 N2=l,3 354 

GRADCO(N2,i »=D(N2, J) 355 

705 DSAVE(N2)=D(N2, J) 356 

RSAVE=0.0 357 

RSUM=0.0$KMAX=0.0 358 

INUM=0 359 

JNUM=0 360 

DO 771 1D=1 ,100 361 

RLAST=RSAVE 362 

RSAVE=RSAVE+D1ST/(10.0J 363 

DU 701 1 = 1 , ICNT 364 

R=0.0 365 

IF ( INUM.EO.IOO) GU TO 701 366 

IF (D( 1 , N .EU.TESr.OR.l .EO. J » GO TO 701 367 

DO 703 Nl=l,2 368 

IF(N1 .EO.l ) DX=DSAVE(N1)-D(N1 , I ) 369 

IF(N1 .E0.2I DY=USAVE(N1 )-D(Nl ,I ) 370 

703 R=R+(DSAVE(N1»-0(N1,1)H'*2 371 

R=S0RT(R) 372 

IF (R.GT.RSAVE) GO TO 701 373 

IF (R.LE.RLASn GO TO 701 374 

IF (ABS( (OSAVEI J)-D(3, I) I ) /R. G T . TANG . QT . A6S ( DS AVEI 3 I -0 ( 3 375 

1 , I n .GT. lOMXI GO TO 701 376 

KADD=-1 377 

IF ( INUM.GE.KMAXl ) 773,7774 378 
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773 00 77i) 110=1 tlNUM 379 

IF (RMAX.ty .RSUIO) ) 77A.775 380 

774 CONTINUE 381 

KAUD=KADD+1 382 

IF IKAUO.GT.O) GO TO 7774 383 

77b CONTINUE 384 

GO TU 772 385 

7774 CONTINUE 386 

IF (R.GT.RMAX) RMAX=R 387 

INUM=INUN+1 388 

RS(INUM)=R 389 

DO 7C7 N4=l,3 390 

DXIVX=AMAXll OXNX.OX) 391 

UYMX=AMAX1 ( DYMX f OY ) 392 

OXi^N=AMiNUOXMN,L)XJ 393 

DYNN=AMIN1 ( DYMN ,0Y > 394 

707 BSAVEIN4, INUK)=D(N4f I ) 395 

701 CONTINUE 396 

771 CONTINUE 397 

772 CONTINUE 398 

IF ( INUN.LT.KMAXl » GO TO 751 399 

DO 712 12=1, I NUM 400 

712 KSUM=RSUM+( 1.0-RS( 12 ) /RMAX I **2*I RS ( 12 » /RMAXJ**S 401 

RWSUM=RSUN 402 

INUM1=1NUM-1 403 

DO 710 11=1 , INUMi , 1 404 

NNUM=I1 405 

JNUM=I 1+1 4SUMS0 1=0.0$ SUMS02=0.0 406 

DO 710 16= JNUM, INUM.l 407 

DO 770 J6=l,3 408 

GKADCCI J6,2 ) = BSAVE( J6, 1 1 ) 409 

GRAUCGI J6,3 ) = BSAVE( J6,'I6) 410 

IF (J6.E0.3I GU TO 770 411 

SUMSai = (GRADCO( J6,2)-GRADC0( J6,1H **2 412 

SUMS02=(GRADCO( J6,3 )-GRAUCO( J6,l )> **2 413 

770 CONTINUE 414 

IF (SUMSai.EO.0.0) GO TO 710 415 

IF ISUMSa2.E0.0.0 1 GO TO 710 416 

RTSMS01 = SClRT(SUMSyi) 417 

RTSMS(.2 = SuRT(SU,MSU2) 418 

SALPHA1=A6S(GRADC0(2,2>/RISMS01 ) 419 

SALPHA2=A8S (GRAOCOI 2,3) /RISMS02) 420 

C ALPHA 1=ABS(GRADC0( 1 , 2 ) / RT SMSQ 1 ) 421 

CALPHA2=A3S (GRADCOI I , 3) /RT SMSQ2 ) 422 

SDIFALP=SALPHA2*CALPHA1-CALPHA2*SALPHA1 423 

IF IABSISDIFALP )-00.1 7365.lt. 0.0) GO TO 710 424 

DG=GRADC0(1 ,1 )*IGRADC0(2,2) *GRADC0(3 ,3)-GRA0C0( 3,2)<‘GRA0C0( 2, 3) ) 425 

1-GRA DCO( 2,1 )*( GRAOCOI 1,2 )^=GRA0CU(3,3)-GR A OCO( 3,2) *GRADCO( 1,3) ) 426 

2+GRADCOI 3, 1 ) *( GRADCO ( 1 , 2 ) *GRAOCO( 2 , 3 )-GR AOCO I 2 , 2 l^GRADCOI 1,3)) 427 

IF (DG.EO.O.OGQ TO 710 428 

GRADIl )=(GRADC0I2,2)*GRADC0(3,3)-GRADC0( 3,2)*GRAJC0I2,3)-GRAD:0I2, 429 

11 ) *(GKADCO( 3,3) -GRADC0I3 ,2 ) ) 4GRADC0I 3, 1 ) * ( 3R AOCOI 2 , 3 )-GRADCO ( 2 , 2 ) ) 430 

2)/DG 431 

GRA0I2 ) = (GRADCO( 1 , 1 ) # ( GR AOCO 1 3 , 3 )- GR ADCO I 3 , 2 ) ) - I GR ADCOI 3 , 3 ) *GR ADCO 432 

11 1 .2) -GRADCU{3 , 2 ) *GRA0C3( 1 ,3 ) ) + ( GR ADCO I 3 , 1 ) ♦ I GR ADC 0 ( 1 , 2 ) -GRADC 3 ( 1 , 433 

23))))/DG 434 

GRAU{3) = (GRADC0( 1 , 1 ) * ( GR AOCQ ( 2 , 2 ) -GR ADCO ( 2 , 3 ) ) -GR AOCO 1 2 , 1 ) * I GR AOCO 435 

11 1.2) -GRADCG( 1 ,3) ) +GRAOCOI 1 , 2 ) *GRADC 0( 2 , 3 ) -GR ADCQ ( 2 , 2 ) »GRADC3( 1,3) 436 
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21 /OG 437 

RW=SQK7 ( (OSAVE( U-BSAVE( 1 ,11 ) )**2*-iDSAVE(2) -BSAVE(2, m )**2)/RMAX 438 

w=SCRT( ( 1 .0-RW )**2*kW##S/RwSUMJ 439 

E( 1 )=W*( 3SAVE( i ,Ii)-CSAVE(lJ) 440 

E(21 = W#(8SAVE(2,m-CSAV£(2l > 441 

E( 3)=W*( BSAVEl 3, I1)-DSAVE( 3) » 442 

00 711 Jl=l,3 443 

BBlJlfU^BBCJl.ll-EUD^CEd )*GRAOI 1 M-E ( 2 » *GRAD ( 2 ) +E ( 3 ) *GRAO( 3 ) ) 444 

1<=(-1.0» 445 

DO 711 J2=l ,3 446 

AA(J1,J2»=AA(J1,J2)+E( J1)*E( J2) 447 

711 CONTINUE 448 

710 CONTINUE 449 

IF (JNUM.EQ.U GO TO 780 450 

CALL MATINV(AA,3.BB,l,OETERM,IPIVUT, INDEX,3, ISCALEI 451 

DO 1700 IA=1,3 452 

IV=LEGVAR(BS( IA,1 n 453 

IF IIV» 751,1700,751 454 

1700 CONTINUE 455 

IF (DETERM.ED.O.O) GO TO 751 456 

DC=83(3,1) 457 

A(J)=BB(1,U 458 

B(JMBa(2,n 459 

IFnuC.EQ.O.OI . AND. { A( J » . EO. 0.0) .AND.( B( J » . EQ.O.On GO TO 751 460 

BB(3, 1)=DC/S0RT(DC**2+A( J ) 2+B I J ) **2) 461 

GU TO 781 462 

780 CONTINUE 463 

A(J)=GRAU(1) 464 

B(J)=GRAD(2) 465 

OC=GRAO(3) 466 

GU TO 751 467 

781 CONTINUE 468 

TXy=(OYMX-UYMN) /< JXMX-OXMN) 469 

B3(3.1)=ABS(BB(3,in 470 

IF (HB(3.1) .LT.CSLQPE) GO TJ 751 471 

A( J )=-A( J )/0C 472 

B( J)=-B( J)/DC 473 

C( J)=-( At J) *DSAVc( 1) +B( J)')‘DSAYE(2)-0SAV£ ( 3) ) 474 

GO TO 700 475 

751 A(J)=TEST 476 

B(J)=TEST 477 

C(J)=TEST 478 

700 CONTINUE 479 

DO 720 1=1, KK 480 

X=KOKE( 1 ) +DIST^FL3AT ( I-l ) 481 

DO 720 J=1,JJ 482 

Y=X0RN(1)-DIST*FLQAT( J-1 ) 483 

RMAX=0.0 484 

RSAVE=0.0 485 

INUM=0 486 

DO 742 L=l,100 487 

RLAST=RSAVE 488 

RSAVE=DIST#FLOAnLI/(10.0) 489 

DO 741 II=1,ICNT 490 

KADD=-1 491 

IF (INUM.GE.KMAX) 776,7777 492 

776 DO 778 13=1, INUM 493 

IF ( iURTI (0(1 , I I )-X)«*2 + (D(2 , II )-Y )#*2) .EO.RMAX) 777,778 494 


23 



777 CONTINUE 495 

KACD=KADO+l 496 

IF (KADD.GT.O) GO TO 7777 497 

778 CONTINUE 498 

GO TO 744 499 

7777 CONTINUE 500 

IF ( A( I I > .£0. TEST. ANU.8( 1 1 I .EJ.TESD GO TO 741 501 

OR = SuRT( (U( 1, 1 1 D(2. I n-Y »**2) 502 

IF IDR.GT.RSAYE) GO TO 741 503 

IF (OK.LE.RLAST) GO TO 741 504 

IF (DR.GT.KMAX) RMAX=DR 505 

INUM=INUM+1 506 

1R(INUM)=II 507 

741 CONTINUE 508 

742 CONTINUE 509 

744 CONTINUE 510 

IF( INUM.EO.OGO TO 720 511 

RD=0.0 512 

DO 730 11=1, INUM 513 

I4=1K(I1) 514 

R=S0kT( CX-U( l,l4))**2 + (Y-0(2,I4))=t=*2)/RMAX 515 

730 RD=KD+ ( 1.0-R) >!'=!=2^‘R<'=!'S 516 

IF(RU.EO.O.O)GU TO 720 517 

DO 731 12=1 , INUM 518 

I4=IK(I2) 519 

R:>=SURT( ( X-D( 1 , 14 ) )**2 + (Y-D( 2. 14) ) ♦S'21 520 

MAT< J, I )=MAT( J , I H-( 1 .0-RS/RM4X)**2*( RS/RMAX) <=*S/RD«' ( C ( I 4) +AI I 4) 521 

16(I4)*Y) 522 

731 CONTINUE 523 

MAT! J , I ) = MAT( J , I ) •■ lADDI 524 

720 CONTINUE 525 

3333 RETURN 526 

END 527 
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Subroutine REGRID 


REGRID first scans the depth matrix generated by GRIDIT to evaluate its maximum 
and minimum values (IMAX and IMIN) as well as the munber of zero (unfilled) elements 
therein. If a border column or row contains all zeros, the matrix dimensions are reduced 
by TESTI. 

Each zero element is filled and checked by REGRID. Further, the slope between 
each pair of adjacent elements in the matrix is checked. If the slope exceeds either 
TANG or IMAX or is less than IMIN, the element is first rejected and then evaluated 
again. It is possible for this routine to reject the entire depth matrix and terminate the 
job. 

SUBROUTINE REGRID 528 

INTEGER AGAIN 529 

DIMENSION 1(60,60) 530 

COMMON /BLKl/ I 531 

COMMON /8LK3/ KM.LM 532 

COMMON /BLK4/ DIST.TANG.M 533 

COMMON /LIMITS/ M I N I , M AXI , I A DD I 534 

600 FORMAT MINIMUM ELEVATION PRIOR TO SMOOTHING = *110) 535 

601 FORMAT (* MAXIMUM ELEVATION PRIOR TO SMOOTHING = *110,/) 536 

602 FORMAT!/,* MINIMUM ELEVATION AFTER I NTERPOL AT I DN ETC = *,110,/) 537 

603 FORMAT {* MAXIMUM EXPECTED SLOPE BET»(EEN T<<0 CONSECUTIVE GRID P3I 538 

INTS = *,F15.5) 539 

604 FORMAT (2016) 540 

605 FORMAT (1015) 541 

607 FORMAT (* UNSMOOTHED DEPTH MATRIX* //) 542 

606 FORMAT (* SMOOTHED DEPTH MATRIX*//) 543 

25 FORMAT (110,110) 544 

6 FORMAT (13,F10.0,F10.5) 545 

5 FORMAT (2014) 546 

IF ( MINI .LT. 1000) 300,302 547 

300 PRINT 301 548 

301 FORMAT!* YOU ARE IN SUBROUTINE REGRID WITH A MINIMUM SCALAR VARIAT 549 

IION WHICH IS TOO SMALL*/* CHECK LAST MINIMUM SCALAR VARIATION IF 550 

2PR0GRAM FAILS TO EXECUTE PROPERLY*) 551 

302 CONTINUE 552 

310 FORMATdHl,* THE FOLLOWING INFORMATION REGARDING YOUR UNSMOOTHED D 553 

lEPTH MATRIX IS FURNISHED*///) 554 

311 FORMAT!/) 555 

312 FORMATdHl,* THE FOLLOWING INFORMATION REGARDING THE SMOOTHED DEPT 556 

IH MATRIX IS FURNISHED*///) 557 

PRINT 310 558 

CALL TESTI 559 

PRINT 311 560 

M=1 561 

MCNT=0 562 

PRINT 607 563 

DO 711 K=1,KM 564 

PRINT 604,( I(K,J) ,J=1,LM) 565 

DO 711 L=1,LM 566 

IF ( I (K,L) .LT.MINI .OR.I (K,L) .GT.MAXI ) I(K,L)=0 567 

IF ( I(K,L).EQ.O) MCNT-=MCNT + 1 568 

IF ( I(K,L).NE.0) 800,801 569 
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800 IF IMCNT.GT.MJ M=MCNT 570 

MCNT=0 571 

801 CONTINUE 572 

711 CONTINUE 573 

PRINT 312 574 

IF (M.EQ.O) PRINT 313 575 

313 FORMAT!* IT HAS BEEN DETERMINED THAT EACH ELEMENT OF YOUR UNSM3QTH 576 

lED DEPTH MATRIX IS FILLED*/* NO ELEMENT WAS FOUND TO EXCEED THE LI 577 

2MITS ESTABLISHED BY YOU*/) 578 

IF (M.NE.O) PRINT 314 579 

314 FORMAT!* IT HAS BEEN DETERMINED THAT EACH ELEMENT OF YOUR UNSM30TH 580 

lED DEPTH MATRIX IS INADEQUATELY FILLED*/) 581 

CALL TESTI 582 

PRINT 311 583 

JCNT=0 584 

ICH£AT=0 585 

THIS SECTION FINOS MIN AND MAX DEPTHS 586 

IMIN=20000 587 

IMAX=-10000 588 

DO 710 K=lfKM 589 

DO 710 L=1,LM 590 

IF !I!K,L).GT.IMAX) 701,702 591 

701 IMAX=1!K,L) 592 

702 IF !I!K,L).LT.IMIN.AND.I!K,L).GT.O) 703,710 593 

703 1MIN=I!K,L) 594 

710 CONTINUE 595 

MAXL=IMAX-1ADDI 596 

MINL=IMIN-I ADD! 597 

PRINT 315 598 

315 FORMAT!* REGRID WILL NOT PERMIT THE INTERNAL ELEMENTS OF I TO EXCE 599 

lEO THE FOLLOWING LIMIT, S*/) 600 

PRINT 600,MINL 501 

PRINT 601,MAXL 602 

700 ICNT=0 603 

AGAIN=500.0 604 

K1=KM-1 605 

L1=LM-1 606 

IF IM.EQ.OI 1001,1 607 

1 CONTINUE 608 

THIS SECTION LINEARLY INTERPOLATES TO FILL MISSING INTERIOR GRID P 609 

DO 20 K=2,K1 610 

DO 20 L=2,L1 611 

IF !I!K,L).EQ.O) 10.11 612 

11 I!K,L)=I!K,U 613 

GO TO 20 614 

10 IF !I!K-1,L).EO.O.OR.I!K+1,L).EQ.O) 12,13 615 

13 I!K,L) = !I!K-l.L)«-l!K+l,L))/2 616 

GO TO 20 617 

12 IF I I!K,L-1).EO.O.OR.I!K,L+1).EO.O) 14,15 618 

15 I!K,L)=! I!K,L-1)+IIK,L*1) )/2 619 

GO TO 20 620 

14 IF !I!K+l,L>.EO.O.OR.I!K,L+i).EQ.O) 8,17 621 

17 I!K,L)=!I!K+l,LH-I!K,L + l))/2 622 

GO TO 20 623 

8 IF ! IIK-1,L) .EQ.O.OR.I!K,L-l).EQ.O) 16,7 524 

7 I!K,L) = !I!K-l,U + I!K,L-l))/2 625 

GO TO 20 626 

16 IF {IIK.L + D.EO.O.OR.IIK-l.D.EQ.O) 18,19 627 
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19 I(«.,L) =( I(K,L<-1> + 1(K-1.L))/2 628 

GO TO 20 629 

18 1F(I(K+1,L».EJ.O.OR.I(K,L-1) .EQ. 0)20,21 630 

21 I(K,L) = ( HK+l.L )+I (K,L-i) )/2 631 

20 CONTINUE 632 

iCNT=ICNT+l 633 

IF (ICNT.EO.M) 23,1 634 

23 CONTINUE 635 

C FILL BLANKS I(1,L) BY LINEAR INTERPOLATION 636 

DO 90 L=2,L1 637 

DO 90 K=l,l 638 

IF (HK,L).EG.O) 91,90 639 

91 IF ( I (K.L4-1 ) .EU.O.OR. I (K,L-1) .EO.O) 90,92 640 

92 nK,L) = (I(K,L+ll + I(K,L-l))/2 641 

90 CONTINUE 642 

C FILL BLANKS 1(1, L) BY EXTRAPOLATION 643 

DO 50 L=1,LM 644 

DO 50 K=l,l 645 

IF ( I (K,L ) .EC.O.AND. ( I (K4-1,L ) .NE.O.AND. I (K*-2,L) .NE.O) ) 51,50 646 

51 J=I(K+1,L)-I(K+2,L) 647 

I (K,L) = I(K + 1,LH-J 648 

50 CONTINUE 649 

C FILL BLANKS I(K,1) BY LINEAR I NTEKPCLAT I ON 650 

DO 30 K=2,K1 651 

DO 30 L=l,l 652 

IF ( I (K.D.EG.O I 31,30 653 

31 IF ( I(K<-1,L).EO.O.OR.IIK-1,L).EO.O) 30,32 654 

32 I(K,L) = ( IIK4-1,L)+I(K-1,L) )/2 655 

30 CONTINUE 656 

C FILL BLANKS I(K,1) BY EXTRAPOLATION 657 

DO 55K=1,KM 658 

DO 55 L=i,l 659 

IF ( I (K,L) . EL.O.ANU. ( I (K,L4-1 ) .Nh.O .AND. I U,L+2) .NE.O) I 56,55 660 

56 J=I(K,L+l)-I(K,L+2) 661 

I (K,L)=I (K,L+1 )+J 662 

55 CONTINUE 663 

C FILL BLANKS llKM,L) BY LINEAR INTERPOLATION 664 

DO 35 K=KM,KM 665 

DO 35 L=2,L1 666 

IF ( I (K ,L) . EG.O) 36,35 667 

36 IF ( I (K.L+1) .EO.O.QR.I ( K,L-1 I .EO.O) 35,37 668 

37 I (K,L)=(UK,L<-1 )+I(K,L-l) )/2 669 

35 CONTINUE 670 

C FILL BLANKS I(KM,L) BY EXTRAPOLATION 671 

00 60 K=KM,KM 672 

DO 60 L=1,LM 673 

IF ( I (K,L) . EG.O.ANU. ( 1 ( K-i ,L ) .NE.O .ANO.I (K-2,L) .NE.O) ) 61,60 674 

61 J=I(K-1,L)-I(K-2,L) 675 

I(K,L)=I(K-1,L)+J 676 

60 CONTINUE 677 

C FILL BLANKS i(K,LN) BY LINEAR INTERPOLATION 678 

DO 670 K=2,K1 679 

DO 670 L=LM,LM 680 

IF I I(K,L) .EG.O) 671,670 631 

671 IF 1 1 ( K + 1 ,L ) .EG.O .OR. I ( K-1 ,L ) . EG.O ) 670,672 682 

672 I(K,L) = (I(K + l,L)4-I(K-i,L))/2 683 

670 CONTINUE 684 

C FILL BLANKS I(K,LM) BY EXTRAPOLATION 685 

DO 65 K=1,KM 686 
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DO 65 L=LM,LM 687 

If ( I (K.LI .EU.O.AND. 1 1 (K,L-1 » .NE.O.AND. It K,L-2» .ME.On 66,65 688 

66 J=1(K,L-1»-I (K,L-2» 689 

I(i<,Ll = I(K,L-l)+J 690 

65 CONTINUE 691 

C AT THIS POINT ALu BLANKS ARE FILLED 692 

1001 CONTINUE 693 

ICNTI=0 694 

IF (JCNT.GT.O) 999,998 695 

998 CONTINUE 696 

DO 997 K=1,KM 697 

DU 997 L=1,LM 698 

IF { nK.LI .GT. INAX .UR. II K, L) .LT. IMIN ) 996,997 699 

996 I(K,L)=0 700 

350 FORMAT!* BE CAUTIOUS OF ELEMENT I ( *I 2* , * I 2* ) * ) 701 

PRINT 350, K,L 702 

ICnTI=ICNTI+1 703 

997 CONTINUE 704 

M=20$JCNT=1 705 

IF (ICNTI.EU.O) GO TO 999 706 

_ GO TO 700 707 

999 CONTINUE 708 

4 FORMAT (2X,2015) 709 

IF (I ( 1 , i ) .GT . I ( i ,2 ) ) 76,77 710 

77 MIN=I!1,1) 711 

GO TO 78 712 

76 MIN=I (1 ,2) 713 

76 CONTINUE 714 

DU 75 K=2,K1 715 

DU 75 L=2,l1 716 

IF ( IIK,L) .GT.MIN) 75,80 717 

80 CUNTINUE 718 

MIN=I(K,L) 719 

75 CONTINUE 720 

MInL=MIN-IADDI 721 

PRINT 602.MINL 722 

If (JCNT.Ey.3) GO TO 305 723 

If (MIN. G£. 1000) 81,303 724 

303 IF (JCNT.EQ.l) 304,305 725 

304 PRINT 306 726 

306 FORMAT!* MINIMUM SCALAR VARIATION IN INTERNAL ELEMENTS OF I MATRIX 727 

1 IS rOU SMALL*/* rilLL TRV TO SMOOTH AND THEN EVALUATE AGAIN*/) 728 

GO TO 81 729 

305 PRINT 307 730 

307 FORMAT!* HAVE TRIED 10 EVALUATE DEPTH MATRIX*/* WILL PRINT I MATRl 731 

IX BEFORE TERMINATING VOUR JOB*/* PLEASE CHECK AND RESUBMIT*//) 732 

DO 320 K=1,KM 733 

PRINT 604,! I!K,U ,L=1,LM) 734 

320 CONTINUE 735 

PRINT 308 736 

308 FURMAT!* POSSIBLE REASONS FOR FAILURE*//* ABSOLUTE VALUE OF MINIMU 737 

IM AND MAXIMUM SCALAR VARIATIONS MAY BE TOO LARGE*/* DISTANCE 3£Tw£ 738 

2EN GRID POINTS NEEDS ADJUSTING*/* INSUFFICIENT CONTROL POINT OENSI 739 

3TY*/* MAXIMUM SLOPE BETWEEN DATA AND/OR GRID POINTS NEEDS ADJUSTl 740 

4NG*) 741 

PRINT 309 742 

309 FORMAT!* PROGRAM EXECUTION TERMINATED WITH A MODE 4 FATAL ERROR*) 743 

AZERO^O.O 744 
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AX^BX/AZERO 745 

A/=AX«AX 746 

82 CONTINUE 747 

PRINT 606 74 a 

DO 9912 K=1 ,KM 749 

PRINT 604,< I(K,U ,L = 1,LM) 750 

9912 CONTINUE 75 I 

REWIND 8 752 

RETURN 753 

81 TAM=TANG 754 

PRINT 603fTAM 755 

PRINT 311 756 

K 2 =KM -2 757 

L2^LM-2 758 

DU 95 K=1.K2 759 

DO 95 L=1,L2 760 

KP1=K+1 761 

KP2=K+2 762 

LPi=L+l 763 

LP2=L+2 764 

T1=ABS( ( l(K,U-KKPl,L))/DiST) 765 

IF ( TAM-Tl ) 100»95,95 766 

100 T2= n (KPi ,L)-I(KP2 .L) J /DIST 767 

I 3-A8S(T2) 768 

IF (TAM-T3) 102,103,103 769 

103 I(KtU = I (KP1,L) I (KPi,L)-I ( KP2,L) ) 7 70 

GO TO 95 771 

102 T4=AB6C ( I( K , LPl I ( K, LP2) ) /DI ST ) 772 

IF (TAM-T4) i05#106,106 773 

106 I(KfU = nK,LPl)+U(KfLPl)-I(K,LP2) ) 774 

GO TO 95 775 

105 KK.D^O 776 

95 CONTINUE 777 

DO 110 K=1,K2 778 

DU 110 L=1,L2 779 

LP2=L+2 780 

KP2=K+2 781 

LPl^L+1 782 

KP 1 =K +1 783 

T1C=ABS( (l(K»t)-I(K,LPl))/OiST) 784 

IF ( I AM-TIO) 200,110,110 785 

200 T20=( I (K,LPn-l (K,LP 2 ) ) /DIST 786 

T30=A6S(T20I 787 

IF «TAM-T30) 202,203,203 788 

203 I ( K,L 1=1 ( K. LPl ) +( 1 { K,LP1)-I ( K,LP2I ) 789 

GO TO 110 790 

202 I(K,L)=0 791 

110 CONTINUE 792 

DO 650 K=1,KM 793 

DO 650 L=1,LM 794 

IF ( I(K,U •GT.IMAX.QR.1(K,L) .LT.IMIN) 651,650 795 

651 I(K,U=0 796 

650 CONTINUE 797 

DO 500 K=1,KM 798 

DO 500 L=1,LM 799 

IF (KK.U.EO.U) 501,300 800 

501 AGAIN=0 801 

500 CUNTINUE 802 

JCNT=JCNT+1 803 

IF lAGAIN.EQ.O) 700,82 804 

END 805 
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Subroutine CKPOINT 


CKPOINT is responsible for averaging the scalar variation of any control point 
which has been observed more than once. 

SUBHLIUTlNh CKPOINT ( X I Ni , Y IN 1 , I UEP I , IPTOK , ICUOE ) 806 

IF (1C00E.EU.9) GO TO 112^ 807 

IF (IPTOK.EQ.O) 103,105 808 

1123 CONTINUE 809 

IUEP2=IDEPSUM/IUtNQM 810 

lDtPSUM=0 311 

IUhNQM=l 812 

GO TO 1133 813 

1113 lOENCM^ IDENGM+1 814 

ID£PSUM=IDEPSUM+IDEP1 815 

IPTOK^-1 816 

GO TO 104 817 

103 IPTOK=-l 818 

IUEN0M=1 819 

X1N2=XIN1 820 

YIN2=YIN1 821 

IDEP2=IDEP1 822 

1DEPSUM=I0EP2 823 

GO TO 104 824 

105 CONTINUE 825 

IF ( YIN1.E0.YIN2. ANU.XIN1.EQ.XIN2) GOTO 1113 826 

1124 IF I IDENOM.GT.l ) GO TO 1123 827 

1133 CONTINUE 828 

TE^'P=XIN2 829 

XIN2=XIN1 830 

X1N1=TEMP 831 

TEMP=YIN2 832 

YIN2=YIN1 833 

YIN1=TEMP 834 

TEMP=IUEP2 835 

IU£P2=I0EP1 836 

IDEP1 = TEP.P 837 

I UEPSUM=lDfcP2 838 

IPT0K=1 839 

IF( IC0UE.E0.9) IPrOK=C 840 

104 RETURN 841 

ENO 842 
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Subroutine TESTI 


TESTI scans the depth matrix for border rows and columns for which the sum is 
zero. If this condition occurs, the appropriate dimensions of the matrix are reduced. 

SUBROUTINE TEST I 843 

REAL MINWRU2.MAXWRD2 ,MINWRD3 ,MAXWRD3 844 

DIMENSION IR0W(2) , IC0L(2) 845 

COMMON /8LK1/ 1(60,60) 846 

COMMON /8LK3/ KM,LM 847 

COMMON /3LK5/ M 1 NWRU2 , M AXWRD2 , M INW RD 3, M AXWR03 848 

COMMON /BLK4/ DIST 849 

101 CONTINUE 850 

DO 100 J=l,2 851 

IR0W(J)=0 852 

DO 100 K=1.LM 853 

IF ( 1( J,K) .GT.O) IR0W( J)=iR0W( J)+l 854 

100 CONTINUE 855 

IF (IR0W(lH*IR0w(2) .GT,0) GO TO 302 856 

DO 300 J^1,KM 857 

00 300 L=1,LM 858 

300 I( J,L)=I( J+1,L) 859 

KM-=KM-1 860 

MAXhRD2=MAXV^K02-DiST 8 61 

GO TO 101 862 

302 J=KM 863 

IR=0 864 

00 301 L=1,LM 865 

IF (i(J,U.GT.O) IK-IR+1 866 

301 CONTINUE 867 

IF CIR.GT.Oi 201,304 868 

304 KM=KM-1 869 

M 1 NWRD2 = M I NyJKD2+D 1ST 8 70 

GO TO 302 871 

201 CONTINUE 872 

DO 200 K=l,2 873 

ICOL(K)=0 874 

DO 200 J=1.KM 875 

IF {1IJ.X).GT.0) ICOL(K)=ICGL(K)+l 876 

200 CONTINUE 877 

IF ( IC0L(l)+iC0L(2) .GT.O) GO TO 402 878 

DU 400 L=1,LM 879 

DO 400 J=1,KM 880 

400 I ( J»L)=I( J,L+i > 881 

LM-=LM-1 882 

M1NWRD3=MINWRD3+0IST 883 

GO TO 201 884 

402 J=LM 885 

1R=0 886 

DO 401 L=1,KM 887 

IF (l(L.J).GT.O) IR=1R4-1 888 

401 CONTINUE 889 

IF (IR.GT.O) 501,404 890 

404 LM=LM-1 891 

MAXWRD3=MAXWRD3-D1ST 892 

GO TG 402 893 
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501 CONTINUE 894 

PRINT 1 .MINWRD2 ,MAXNR02 ,MIN<JRD3fMAXWRD3 895 

1 FORMATI/, * 896 

1 ADJUSTED BOUNDARIES UF DEPTH MATRIX*//* MINIMUM Y*E16.8* MAXIMUM Y 897 

1*E16.8,/,* MINIMUM X*E16.8* MAXIMUM X*E16.8) 898 

RETURN 899 

END 900 
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Subroutine CONVERT 

Program variables are scaled to plot dimensions and converted to inches by 
CONVERT. 

SUBKUOTI.Nfc CuNVEKT(X) 901 

CGMMCN /3LK2/ K JR E ( 4 ) , K0R(V { 4 ) , i F YP , S CL E 902 

CUMMQnI/ 3LR14/ UNIT.CNVTJIN 903 

X=X*uNI r/( SCLE^CiMVTUIiM) 904 

RETURN ^05 

end 906 
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Program CKZFIT 


CKZFIT is the second of three overlays which are subordinate to COUPLE. 

CKZFIT establishes the X- and Y-coordinates of each element within the depth matrix 
and stores each within the appropriate array table and reorganizes the I matrix into a 
one-dimensional array. Each of the array tables and each control-point X- and 
Y-coordinate are part of the argument listed for the CALL statement to DISCOT, a 
system routine. DISCOT uses linear interpolation to find the scalar variation predicted 
by the depth matrix at each. X and Y control-point position. The real and predicted 
control-point scalar variations are stored in turn by CKZFIT. Subsequent to reading all 
the control points, subroutine SUBl is called. 

Ov/fcKLAy (LINK., 5,3) 907 

PKUuKAM CKIFIT 908 

01 Mdi'Ji.ION XT Ad ( o J ) , ZTA3 ( 6C ) , riAB( 3b0 0) , I ( oO ,60 ) ,ARG( 1005,2 ) 909 

OIMtNSlCN K.U.Rt (4 ) ,K.UR!\i ( A) 910 

DIMENSION YYT AO (a 5,60) 911 

EuUlVALEN'Cfc IYYTAb,YTAB) 912 

CCMMCN /3LK1/ I 913 

LUMMbi', /BLK2/ KORcKORN 914 

CCMMLN /faLK3/ NM,LM 915 

CCMMUN /OLKA/ 0 L ST 916 

CCMMQ.N /LIMITS/ M I N I , M Ax I , I A DU I 917 

Cu.MMON /LINKUC/ LINK, RECALL 918 

PRINT 1 919 

i KJRMAT ( i.Hl , =!' YdUR SMUuTMEu DEPTH MATRIX, J, IS PRINTED ABOVE.*/* A 920 

IT IniS POINT THE 6CALAR VARIATION PREUlCTcD BY i AT THE LOCATION 921 

2UE EACi CuNTK'OL POINT*/'* IS COMPARED WITH THE SCALAR MAGNITJDc OF 922 

3THt Control point.*/* a plot for FREDUENCY QF occurrence vs DISCRE 923 

4PANCY FOLLOWS AS PART OF THIS LISTING*//////) 924 

REWIND 7 925 

YMAX^FLOAT { KORN ( 1 ) ) -FLOAT ( KM-1 ) *Di ST 926 

XM1N=FL0AT( KuKEI i ) ) 927 

1C,\I = 0 928 

I l=C 929 

Ju) luO L = 1,CM 930 

XTAb(L)=XFlN + Ol,ST*(FLJAr(L-l ) ) 931 

i 1 = ( 60) *( L-1 ) 932 

I2=11+KM+1 933 

JU ICO 0=1, KM 934 

IF (L.bO.l) ifAtJ( J) =Y,MAX + 0IST*FL0AT( J-1 ) 935 

Ii=12-J 936 

YTABI 11 )=I ( J,L) 937 

100 CONTINUE 938 

LM1=LM*1 939 

KMi=KM+l 940 

Du 101 N=Khl,60 941 

ZTAtt(N) =ZTAB( K.M) +01 ST* ( N-KM) 942 

OU 101 J=1.LR 943 

K= ( J- 1 ) * I 60 ) +KM 944 

101 YYIAu(N,J)=YTaD(K) 945 

DO 102 N = l,KPi 946 

K=(LM-i )*(6C)+N 947 

DO 102 J=LM1,60 948 

XTAB( J) =XTAB( LM) +DIST* { J-LM) 949 
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102 YyTAd(N, J)^VTAl5(K) 950 

200 CONTINUE 951 

CALL RECIN( 7,1 ,KKK, iCUDh, I X, lY, IZt IDISn 952 

IF (hOF,7> 201,202 953 

202 IF( ICQDE.GT .0) GJ TO 201 954 

XA=IX 955 

ZA=IY 956 

CALL 01 SCOT { ZA,XA, ZTAB, YTAB, XTAG,111 , 360 0 , 60 , AN S I 9 57 

ANS=ANS-FLUAT( I AOU I ) 958 

ICNT=ICNT+1 959 

ARGI ICNT,l)=IZ/{ 1000.0) 960 

ARGt ICNT,2)=ANS/( 1000.0) 961 

GO TO 2C0 962 

201 CONTINUE 963 

CALL SU31 ( ARG, ICN J ) 964 

RETURN 965 

ENU 966 
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Subroutine SUBl 


SUBl computes the discrepancy between the predicted and true scalar variations. 
Each discrepancy is then categorized according to its magnitude as described under the 
heading "Numerical Techniques" and is counted along with those discrepancies which 
fall within the same category and which precede it in storage. SUBl, after processing 
each stored variable, plots a frequency distribution curve of the discrepancies as part 
of the listing. The plot is scaled to page size and begins with the first category for 
which the ordinate is not zero and continues until all the nonzero ordinates have been 
plotted. The total number of discrepancies which are found to fall within each category 
is printed to the right of each ordinate in the event the user may be interested in its 
magnitude. Further, the standard deviation and the mean value of all the discrepancies 
are determined. 


SUbKUUTINfc SUbUARG. NOPT) 967 

DIMfcNSlQN ARGiluOO,2),CLASS2(201),DIFF(1000) 968 

UEMENSIO.NJ CLASS(20U ,C0JNT(2Jl),PERCcNT(231) rIM(2) 969 

OIMtASIUN lORD ( li ) . I ABSC( 201 ) ,L JRO ( 201 ) , A3SGJ 201) , XM(2I,YM1(3) 970 

I TAPfc=5LTAPt7 971 

CNTMX=-200.0 972 

REWIND 7 973 

SUM=SUMSy-=81 = 62 = 0.D 974 

ISUM^O 975 

UU 20 1=1, 2C0 976 

20 CLASS ( i ) = ( I-lOl )-i^( 0.02 ) tC. 01 977 

OU 40 1=1,201 978 

40 COUNT ( I )=0 979 

00 130 J=1,N0PT 980 

IF (ARG( J,1 ) .EO.IOOOCOO.O. OR. ARG(J ,2). E3. 1000000.0) GO TOISO 981 

UiFFER=ARG( J,2)-ARG( J,l) 982 

SUM=SUM+U1FFER 983 

ISUM=ISUM+1 984 

SU)^SC = SUMSQ + UI FFER**2 985 

N=J+1 986 

IF ( J.E3.NCPT)GJ TO 45 987 

UO 30 I=N,NGPT 988 

IF (ARGil ,1 ) .EO.IuClCOO.O.QR.ARGII ,2).E0.10000D0.3) GO TO 30 989 

0AKG1=ARG( I , 1 )-ARG( J, 1 ) 990 

0ARG2=ARG( I ,2) -ARG( J ,2) 991 

0ELAkb=CAkG2-0ARi^l 992 

30 CUMINUE 993 

130 CONTINUE 994 

45 CONTINUE 995 

AV=SUK/FL0AT( I SUM) 996 

SI GMAS0 = SUMS0/FLU AT ( ISUM )-AV^i=*2 997 

DO 75 1=1, 2u0 998 

CUUNT(I)=0 999 

75 CLASS2I I ) = ( i-101)‘*'(0.u2)+-C.01 1000 

DO 131 J=1,NGPT 1001 

IF I ARij( J , 1 ) .bO.lOuOCOO. O.OR .ARGI J ,2 ).£0. 1 000000.0) GO T0131 1002 

DIFFER = ARG( J ,2 ) -ARGI J , 1 ) 1003 

81=81+1 DIFFER-AV)**3 1004 

82 = B2+(UIFFER-AV)<‘*4 1005 
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Ulf^H JJ=CIFFEi<-AV 1006 

DO 77 1=1,199 1007 

IM DIFF{ J ) .Gfc ,CLASS2( n , AND. DIFF( J » ,LE ,CLASS2( I +U )G0 TJ 78 1008 

77 CUMINUE 1009 

GO ID 79 1010 

78 COUNT! H-1 ) = Cau;>lT( 1 + 1 )+l ,C 1011 

79 CONTINUE 1012 

131 CONTINUE 1013 

DU 80 1=1, 2JC 1014 

IF (COUNT! 1 ) .GT. CNTMX) CNT '1X=C0UNT ! I) 1015 

80 CLASS2!1)=(CLASS2( 1)+0.01)*1000.0 1016 

lSCALe=!FlX!CNTMX)/10C+l 1017 

JLAST=0$ JFAX=0 1018 

DU 70 1=1,200 1019 

IA8SC ! I )=CLASS2 ! I ) 1020 

LOKD! I) = 1FIX(CUUNT(I)/FL0AT( 1 SCALE) ) 1021 

IF (LORD! I ) .GT .0) JMAX=JMAX + 1 1022 

IF (l.LE.ll) IURD(1)=1-1 1023 

lF{l.LT»96,OR,l.GT.106)AbSC(I)=2H 1024 

70 CONTINUE 1025 

ABSC! 96)=2H D 1026 

A3SC! 97)=2ri 1 1027 

ABSC(98 )=2H S 1028 

ABSCI99 )=2ri C 1029 

ABSC(100)=2h R 1030 

AUSC(101)=2ri E 1031 

ABSC! 102 )=2H P 1032 

ABSC! 103)=2H A 1033 

ABSC! 104)=2H N 1034 

ABSC! 105)=2ri C 1035 

ABSC! 10b)=2H Y 1036 

1START=10H 1037 

IEND=1SCALE<'10 1038 

PRINT 301 , ( lORU ( 1 ) , 1 = 1 , 1 U , lEND 1039 

301 FORMAT(52X,*FKEOUENCY UF LCCUKR E NC E* / , 3 X , 1 1 I 1 0 , *X >1= I 4* MUMBER«=) 1C40 

ISYMB=1H* 1041 

1SYNU1=1H 1042 

DU 20C I=i,2C0 1043 

J=LOKU( 1 ) 1044 

J0=1F IX (COUNI ( I ) ) 1045 

IF ( J .GT.O.ANU.jLASr ,Ed,0) JLAST = 1 1046 

IF (J.EU.O) GU TO 303 1047 

JMAX=JMAX-1 1048 

JJ=J+1 1049 

PR INT 300, ABSC ( I ) , 1 ABSC ( I ) , ( I SY,'3B,i< = i ,J ) , ( 1 SYMBl , JK= J J , 1 30 ) 1050 

300 FORNAT! 1H+, A2, 110, 112A1 ) 1051 

GU TO 30A 1052 

303 CONTINUE 1053 

IF (BLAST. EU.O. OR. UMAX. LE.C) GU TO 200 1054 

PRINT. 300, ASSCd ) ,IAKSC( 1) 1055 

304 CONTINUE 1056 

PRINT 302,00 1057 

302 F0RMATI115X, ilul 1058 

200 CUNT 1 NUE 1059 

PRlNT2,AV,SlGi'3ASO 1060 
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2 FURNAn* ALL CLNTROL POINTS HAVE BEEN CUMPAKEJ AND THE RESULTS ARE 1061 

1 SHUwN*/* IT NAS UETLRNINEO FROM THE C JlNP AR I SON S THAT THERE NAS A* 1062 

2/* MEAN LIFFERENCt =*E16.8,/<' AND ST DEV =*E16.8, /,♦ THE DE 1063 

3PTH MATRIX rtlLL NOW HE PLuTTED*/=i= A TWO DIMENSIUNAL CONTOUR CHART 1064 

4WILL BE PREPARED NlIH A SUITABLE BRIO*/* EACH SRID LINE REPRESENTS 1065 

5 AN Integral multiple of DISPBHDISEE input instructions > */* and T 1066 

6HE CLNTGURS ARE POSITIUNEU RELATIVE TO THIS GRID*/* PLOTTING BEGIN 1067 

7S WITH THE UPPER LEFT HAND CORNER OF THE DEPTH MATRIX*/) 1068 

RETURN 1069 

END 1070 
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Program AUTOCON 


AUTOCON is a very complex program which contains the necessary decision- 
making logic for managing the indexii^, the computation, and the plotting of the contour 
line segments which cross each pair of adjacent triangles obtained from each 3 by 3 sub- 
matrix contained in the depth matrix. The basic concept is described in the section 
"Contour Plots." Further, AUTOCON establishes first a grid of equally spaced lines, 
each of which is separated by the scaled equivalent of DISPGRD (an input variable) , and 
then it positions the plotted contours within the grid relative to the upper left-hand 
corner of the grid. Each grid line represents an integral multiple of DISPGRD. Thus, 
X-Y coordinates of the upper left-hand corner of the grid are the truncated values of 

X = Minimum X/DISPGRD 


and 


Y = Maximum Y/DISPGRD 


AUTOCON plots first the grid and establishes the origin of the plot, after which the 
first two adjacent triangles in figure 2 are processed and then plotted. The program 
indices are then reset for the processing and plotting required for the next two triangles 
in figure 2. This process is repeated until the 3 by 3 submatrix is completely plotted, 
after which the indices are reinitiated and a new 3 by 3 submatrix is considered. 

AUTOCON possesses one option which permits the depth matrix to be plotted at a 
finer grid spacing. If SPP (an input variable) is greater than one, the 3 by 3 submatrix 
is subdivided so that a 2SPP + 1 by 2SPP + 1 submatrix is processed. The indexing is 
automatic and SPP cannot exceed 8. 


OVERLAY (LINK, 3, 2) 1071 
PROGRAM AUTG CON 1072 
REAL MINWRD2,MAXWR02,M INWR03,MAXRP.D3 1073 
DIMENSION IDEPH(160,3) 1074 
DIMENSION DEPTH(17,17) 1075 
DIMENSION IEL(60,60I 1076 
COMMON /BLKl/ lEL 1077 
COMMON /QLK3/ MAT, MIT 1078 
COMMON /8LK5/ M INWRD 2, M AXWRD2 ,M INW RD3 ,MAXWRD3 1079 
COMMON /8LK6/ SP, SPP , LAST ,AP , LINEUNT 1080 
COMMON /BLK7/ N, NN, P XI , PYl , XMAX ,YMAX, XMAY, YMAY , C I 1081 
COMMON /BLK8/ SQX,SQY 1082 
COMMON /BLK9/ PCS1X( 20) ,P0S1Y( 20) ,CCNT (20) 1083 
COMMON /BLK10/P0S2X( 20) ,P0S2Y(20) ,C0NTA( 20) 1084 
COMMON /BLK11/P0S3X(20) ,P0S3Y(20) ,C0NT6(20) 1085 
COMMON /BLK12/P0S4X( 20) ,P0S4Y(20) ,C0NTC(20) 1086 
COMMON /8LK13/PGS5X(20) ,POS5Y(2O),C0NTO(2O) 1087 
COMMON /LINKDC/ LINK, RECALL 1088 
CALL CALCQMP 1089 


39 



CALL LEROY 1090 

MINWR02 = MINwFD2/FL:)AT(L INEUNT) 1091 

MIiMWRD3=MINhK03/FLJAT(LINEUNT) 1092 

MAXWR02=MAXWR02/FL3AT(L INEUNT) 1093 

MAX WRD3=MAXWR03/FL0A TIL INEUNT) 1094 

00 9913 I=lt3 1095 

00 9913 J=l,160 1096 

9913 10EPH(J,I)=0 1097 

DO 9914 1=1,17 1098 

00 9914 J=l,17 1099 

9914 DEPTH! I,J)=C.O 1100 

1000 CONTINUE 1101 

Y20=(MAXV«R02-FL0AT( I FI X(MAXWRD2 ) ) ) 1102 

X20=(MINWR03-FL0AT( I F I X I M IN WRD3 ) > ) 1103 

IF (Y20.0T.0.0) Y23=l.C-Y20 1104 

IF IY20.LT.0.0) Y20=ABS(Y20) 1105 

IF (X20.lt. 0.0) X20=X20+1 1106 

Y20=Y20*AP 1107 

X20=X20-AP 1108 

ILINE=IFIX(NiAXWRD3)- IF I X ( M I NWRD 3 ) 1109 

JLINE=IFIX(MAX»JRD2)-IFIX(MINWRD2) 1110 

LINE CNT=0 1111 

258 IF(Y20+{ MAXWR02-MINWRD2-FL0AT{ JLINE) )*AP+0.5.GE.0.0) JLINE=JLINE+ 1112 
11 1113 

IF (Y20+(MAXhPD3-MINWRD3-FL3AT( ILINE) )4=AP+0.5.GE.0.0) ILINE=ILINE+ 1114 
11 1115 

LINECNT=LINECNT+1 1116 

IF (LINECNT.lt. 2) GO TO 258 1117 

X10=0.0 1118 

IF (X20-0.4.LE.0.0) X10=X10-AP 1119 

IF (X20-0.4.LE.0.0) IL I(mE=I LINE«-1 1120 

Y10=-FLnAT(JLINE)<=AP 1121 

CALL GRID (XlOtYlO, AP,AP,ILINE, JLINE) 1122 

UV=1.0/AP 1123 

OkIGlN=FLOAT( I F I X ( M I NWRD3 ) ) 1124 

IF (X20-0.4.LE.0.0) ORI GIN=ORIGIN-DV*AP 1125 

XDIST=ILINE*AP 1126 

yOIST=JLINE*AP 1127 

CALL AXtS(X10,Y10,0.0,XDI8T, ORIGIN, DV,AP,DV,1H ,0.1,-1) 1128 

URIGIN=FLOAT( I F I X (M AXWK02 ) ) -FLO AT ( JL I NE ) 1129 

CALL AXES(X10,Y10, 90. ,YDIST, ORIGIN, DV,AP,DV,IH ,0.1,+1) 1130 

Y2^=V20*(-1. J) 1131 

CALL CALPLT (X20,Y2C,-3) 1132 

9942 CONTINUE 1133 

XMAT=FLJAT( (Nll-1 )/2)*SP*2 1134 

XMAX=XMAI 1135 

XMAY=-XMAX 1136 

YMAX=FLf)An (MAI-1 )/? )=«=SP*(2.J ) 1137 

YMAX=>'MAX + C.05 1138 

YMAY=-YMAX 1139 

PY1=0.0 1140 

PX1=C.0 1141 

SP=SP/S?P 1142 

SJX=SP 1143 

S(jY=-SO 1144 

SPX=SP 1145 

SPY=-SP 1146 

IF (MT.LT.2C) N0C=1 1147 
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IF (Hir .£Q. (MIT/20)*20) NOC=Mir/20 1148 

IF (Mir.GI. <M1T/20)*20> NGC = MIT/23 + l U49 

21 K=0 115.) 

KK=0 1151 

NX=NuC*20 1152 

1S=SPP 1153 

850 K=K+1 1154 

IF (KK) 8850.8851,8850 1155 

8850 JJ=K-1 1156 

JK=K+1 1157 

DO 8852 1 = 1 , MIT 1158 

lUbPHI I ,JJ ) = IDcPH{ I , JK) 1159 

8852 CONTIMUE 1160 

8851 CO,\iTl,NUE ll6l 

779 CONTINUE 1162 

IF (KK-MAT) 9779,612,812 1163 

9779 CONTINUE 1164 

KKK=KK+1 1165 

DO 9933 J=1,MIT 1166 

IUEPHU,IO = IEL(KKK, J) 1167 

9933 CONTINUE 1168 

7bO KN=KK+1 1169 

782 IF (K-3)9850,784, 784 1170 

9850 K=K+1 1171 

GO TO 8851 1172 

784 J=1 1173 

K=1 1174 

698 IX=1 1175 

IY=1 1176 

L=1 1177 

LL=L+IS 1176 

LLL=LL+IS 1179 

M=1 1180 

MM=M+IS 1181 

llb2 

*^=1 1183 

699 0EPTH(L,M)=1UEPH( J,K) 1184 

700 IF ( IDhPH ( J+1 ,(< I ) 3 ,4 ,3 1185 

3 IF ( I Dt P i I ( J , K ) ) 5 , 4 , 5 1 1 8e> 

4 IF ( J .LE .MI T .AND. lUtPHU , K. ) . NE .0 ) GO TO 5 1187 

SUX=SUX+( (2.0) 4SPP‘«=SP) 1188 

IF(J-1)5,6,7 1189 

7 JCKl=J/2 1190 

JCK2=JCt<il*2 1191 

IF (J-JC;<2) 8,9,8 1192 

8 J=J+2 1193 

GO 10 10 1194 

9 J=J+1 1195 

GO TO 10 1196 

6 J=3 1197 

10 JC=JC+2 1198 

IF (JC-NX+1) 598,1851,1851 1199 

1851 SOY=S JY-({2.0)*SPP*SP) 1200 

K=1 1201 

SOX=SPX 1202 

60 TO 1852 1203 

5 XlNCl=iDEPri( 0+1 ,K)-ICEPH( J,K) 1204 

XINC1=XINC1/SPP 1205 

XlNC=AdS(XINCl ) 1206 
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703 L=L+1 1207 

11- (L-LL) 704,703,706 1208 

7C4 IF (XI, NCI) 701,701,702 1209 

701 DEPIHIL ,N)=DEPTH(L-1,M)-XINC 1210 

GO TO 703 1211 

702 ULPTri(L,M)=DtPTH(L-l,M)+XINC 1212 

GO TO 703 1213 

7Gb J=J+1 1214 

LL=LLL 1215 

IF (L-LLL) 699,699,706 1216 

706 JC=J 1217 

J=J-2 1218 

L=1 1219 

K=X+1 1220 

K=M+IS 1221 

IF (M-KMM) 707,707,708 1222 

707 LL = L + I:> 1223 

GO TO 699 1224 

708 KC=K-1 1225 

L=i 1226 

M=1 1227 

LL=L+IS 1228 

798 yiNCl=J£PrH{L,NM)-L)LPTi-l(L,M) 1229 

YIN01=YINC1/SPP 1230 

YlNC.= AdS(YlNCl) 1231 

1803 M=M+1 1232 

IF ((--M.’l) 1 804,1805,1806 1233 

1804 IF (YlNCl) 1801 ,1801 , 1802 1234 

1801 DLPTH(L,M)=0E9Iri(L,M-l)-YINC 1235 

GO TO 1803 1236 

1802 lJEPIh(L,M) = DcPTH(L,M-l). + YINC 1237 

GO TO 1803 1238 

1805 L=L+1 1239 

M=M-IS 1240 

IF (L-LLL) 798,798,799 1241 

799 M=MM 1242 

1243 

L=1 1244 

GO TO 798 1245 

1806 J=JC 1246 

IN!0LX=1 1247 

M2=2 1248 

L2=2 1249 

Ll=2 1250 

Ki-2 1251 

Jl=l 1252 

J2=l 1253 

K2=l 1254 

Ml=l 1255 

100 CONTINUE 1256 

00 9912 1=1,20 1257 

CUNT ( I )=0.0 1258 

CONTA(I)=0.0 1259 

CQNT8(1)=0.0 1260 

CUNTC(I)=0.0 1261 

CONTum=0.0 1262 

POS1X(I)=O.C 1263 

PUS2X(I)=C.O 1264 
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PUS3xm=0.C 1265 

POS4X(1»=0.0 1266 

PUS5X{I)=0.0 1267 

POSlVm^O.O 1268 

POS2Y<I)=0.0 1269 

PU63Y(I)=0.C 1270 

POS4Ym=Ci.O 1271 

PDi.5Y( I )=0.0 1272 

9912 CUNTINUE 1273 

1=1 1274 

11=1 1275 

111=1 1276 

14=1 1277 

15=1 1278 

CALL LlNEl (OEPTHILI tL2) tOEPTH(Ml,M2 ),CI , IMDEX, 1 ,01) 1279 

CALL LIIME2 (DEPTH) Ml , M2 ) .DEPTH (Ji,J2 ), Cl , INDEX. II . D2 ) 1280 

CALL L1NE3 ( OEP TH ( L 1 ,L2 ). UEP TH { J 1 , J 2 I . C 1 , I NDE X . 1 1 I , 03 ) 1281 

CALL LINE4 ( DE P TH ( J 1 , J2 ), DEP TH ( K1 . K2 ). C I , 1 NDE X , 1 4 .04) 1282 

CALLL1NE5 ( DEP T H ( L 1 , L 2 ), OE P TH I M . i^2 ), C I . 1 NDE X , 1 5 .05) 1283 

IF (01) 101,103.103 1284 

101 IF (D2) 133,114,114 1285 

114 N=1 1286 

NN=1 1287 

134 IF (CONT(N)-CONTA(NN) ) 131,132,133 1288 

131 N=N+1 1289 

IF ( CENT ( N ) ) 1 34, i 33, 134 1290 

132 IF (CJNT(N) )1132, 133, 1132 1291 

1132 CALL PLOT ( PU 32 X ( NN ) . P'J S 2 Y ( N N ) , PCS 1 X (N) , P OS 1 Y ( N ) , ^ 0S2 X ( NN ) , 1292 

1 PGS2 Y(NN) ,CUNTA(NIM) , + l ) 1293 

IF (N) 133.133,1134 1294 

1134 IF (CLNTA(KN) ) 13‘t, 133,1 34 1295 

133 IF (03) 338,140, 140 1296 

140 N=1 1297 

NN=1 1298 

153 IF ( CONTB (N )-CUNTA( NN ) ) 149 , 1 30, 138 1299 

149 IF (CCNTB(N)) 152,138,152 1300 

152 NiN=NN + l 1301 

IF (CONTA(NN)) 153,138,153 1302 

150 IF (CQNT8(N)) 1130,138,1150 1303 

1150 CALL PLOT ( P0S3X ( N ) ,PUS3Y (i'J ) ,P JS2X (NN) , P032Y (NN) , P0S2XI n)N ) , 1304 

1 P0S2Y(NN) .CUNT A(Ni'j) ,-l ) 1305 

GO TG 153 1306 

338 N=1 13(j7 

NN=1 1308 

350 IF (CONTb(N)-CUNrA(NN) ) 349,348,349 1309 

349 N=N+l 1310 

IF (CONTB(N)) 350,538,35C 1311 

348 IF (CUNTB(N)) 1348,538,1348 1312 

1348 CALL PLOT ( P0S2 X ( NN ) , PQS 2 Y ( NN ) , PUS 3X (N ),PJS3Y(N ),PUS2X(NN), 1313 

1 PUS2Y(NN) .CONTA(NN) ,-l ) 1314 

IF (CQNTB(N)) 350,538,350 1315 

538 N=1 1316 

NN=1 1317 

542 IF (CGNT(N)-CGNTB(NN))138,540,138 1318 

540 IF (CONT(N)) 539,138,539 1319 

539 IF (CCNTB(NN)) 541,133,541 1320 

541 CALL CALPLT ( PuS I X ( N ) , PUS 1 Y ( >1 ) , 3 ) 1321 

CALL CALPLT ( P0S3 X ( NN ) , PCS 3 Y ( NN ) , 2 ) 1322 

N=N4-1 1323 

NN=NNH 1324 
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60 TO b^tZ 1325 

103 IF (C2» 162tldlil81 1326 

162 N=1 1327 

NN=i 1328 

184 If (CGMT (M-CJ.n:TA(NN) ) 181.182,188 1329 

183 M=N+1 1330 

If (CONT(N)) 184,181,184 1331 

182 IF (CONI(N)) 1182 ,181 ,1182 1332 

1182 CALL PLOT ( PuS 1 X ( N I , PuS 1 7 I N I ,P082X ( N N) , PQS 2Y ( N.M I , P3S2X ( NN ) , 13 33 

1 PUS27(NN) ,CONIA(NN) ,+l) 1334 

If (IM) 181,181,1184 1335 

1184 IF (CL, \'TA(N.M)I 184,181,184 1336 

181 If (031 18/, 438, 438 1337 

187 N=1 1338 

h.\=l 1339 

200 If (CCNraiiMI-CQMAlNM 1 138,198, 199 1340 

199 NN=.'lf\l + I 1341 

IF (CC:\iTA(N(M 1 200,138,200 1342 

198 If (CLnT3(M) 1 198,138,1198 1343 

1198 CALL PLOT ( PCS 3 X ( N 1 , PUS 3 Y ( M 1 , PUS2 X ( N N ) , P 0 S2 Y ( \ M 1 , POS 2X ( MN I , 1344 

1 PCS 2 Y(U.\) ,CuUIA(NN) ,-n 1345 

If (CGNTOINI) 200,138,230 1346 

433 N=1 1347 

CN=1 1348 

450 If (CONTB(N)-CO,\TA(NNl ) 538,4^,8,449 1349 

449 ,\i=.N+i 1350 

IF (CLMlBlMi) 45j,538,45C 1351 

443 If (CCimT3(.\I) ) 1 44 8 , 538 , 1448 1352 

144a CALL Plot ( PCS2X(UN» ,PJS2Y (\ 4) ,PU8 3X (M ),PJ83Y(M » , P US2 X ( N:M ) , 1353 

1 PUS2Y(,‘'i\) ,CC,\tA('4INi) ,-l) 1354 

IF (CLNTB(N))460,a38,45C 1355 

138 If (031 201.203,203 1356 

201 IF (041 227.200,206 1357 

206 N=1 1358 

NN=1 1359 

226 If (CC,\T8(N )-C0NTC( MiSi) ) 223,224,227 1360 

223 N=N+1 1361 

If (CCNTalUn 226,227,226 1362 

224 IF (COiMTdlNII 122 h ,?? 7 , 1224 1363 

1224 CALL PLOT ( PUS 4X { u,\ ) , f US4 Y ( N J ) , PU6 3X ( vi I , POS3 Y( 'J ) , P J S4X ( 9N I , 1364 

1 PLS4Y(,\M .Cu.OTC(iMN) , + l ) 1365 

IT (N) 227,227,222b 1366 

227 IF (U5) 228,230,230 1367 

230 N=1 1368 

NN=1 1369 

242 IF (CCNT0(N>-CJNTC(NNI ) 2o9, 240,241 1370 

239 N,\=NC + i 1371 

IF (CGNTC(NM)) 2h2, 248, 2^*2 1372 

240 It- (CL4T!0(.4)) 1^40,243, 1240 1373 

12+J CALL PLOT (PUS5X(N) ,P0S5Y(M ,P0S4X(NN) ,PJS4Y(NN),PUS4X(NN) , 1374 

1 PUS4Y(N\) ,CUMTC(NM) ,-l ) 1375 

II (CLiMTCd'Nn 2^,2 . 243, 242 1376 

228 N=1 1377 

NC=1 1378 

25a If (CUMTU(N)-CUNrC(:\iM I 1252,253,241 1379 

252 N=N+1 1380 

IF (CCUT0(Nn 255,241,255 1381 

253 IF (CCiNTUINI ) 125 3,241,1253 1382 
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1253 CALL PLOT { PO S 5X i N ) , PUSS T ( N ) , P0S4X ( NN) , P OS^T ( NN ) , P0S4X ( ,'IN ) , 1383 

1 POS4y<NN) .CONTCdMN) ,-l ) 1384 

IF (CONTC(NN)) 255,241,255 1385 

1252 IF ICCNTDIN)) 252,241,252 138fa 

241 N=1 1387 

NN=1 1388 

543 IF (COIMTBINI-CONTUINIM) } 243,54^,243 1389 

544 IF (CJNI8(N}> 545,243,545 1390 

545 IF ICONTUINM) 546,243,546 1391 

546 CALL CALPLT ( P0S3 X ( N ) , P0S3 Y ( N) , 3 ) 1392 

CALL CALPLT (P0S5 X ( NN ) , PDS5 Y ( NN ) , 2 » 1393 

N=N+1 1394 

NN=KN+1 1395 

GO TO 543 1396 

203 IF (04) 256,275,275 1397 

256 N=1 1398 

NN=l 1399 

1274 IF (CUNTB(Nl-CONTC(NN) ) 275,276,277 1400 

277 N=N+1 1401 

IF (CCiMTBIN)) 1274,2 75,1274 1402 

276 IF (CCNr8((M)J 12 76, 275, 1276 1403 

1276 CALL PLOT ( PUS4 X ( IMN ) , P JS4 Y ( N J , POS 3X (N| ) , PUS 3 Y ( 'J ) , P0S4X ( NN ) , 1404 

1 P0S4Y(NN) ,CONTC(r>jiM) ,+l ) 1405 

IF (M 275,275,1274 1406 

275 CONTINUE 1407 

IF (D5) 278,260,280 1408 

278 N=1 1409 

NN=1 1410 

291 IF (CUNTOtN)-CaNrC(NN ) ) 243,289,290 1411 

290 KN=NN+1 1412 

IF (CQNTC(NN)I 291,243,291 1413 

289 IF (CCNTO(N)) 1 289 , 2 43 , 12 89 1414 

128? CALL FOOT ( PUSS X ( N ) , P0S5 Y ( N ) , PUS4X ( N N) , P0S4 Y ( NN ) , PQS4X ( NN ) , 1415 

1 PUS4Y(NN) ,CUNTC(NN) ,-l ) 1416 

IF (CCNTC(NN)) 291,243,291 1417 

280 N=1 1418 

NN=1 1419 

302 IF (CUNTD(N)-CO(MTC(NN) ) 241, 300,301 1420 

301 N=N+1 1421 

IF (CCNTD(N)) 302,241,302 1422 

300 IF (CGNTO(N)) 1 3 'JO ,24 1 , 1 3 00 1423 

1300 CALL PLUT ( PUS 5 X ( N ) , PUS5 Y ( N ) , PUS4X ( N N) , PUS4 Y ( N N > , P J S4X ( NN ) , 1424 

1 PUS4YCMN) ,CCNTC(NN) ,-l ) 1425 

IF (CCNTD(NN)) 302,241,302 1426 

243 INDEX= INOEX+1 1427 

IF (INDfcX-4) 800,800,801 1428 

800 GO TO ( 604,805,806,807) ,I'(UtX 1429 

804 Jl^Jl+2 1430 

J2=J2-2 1431 

K1=K1F3 1432 

K2=K2-1 1433 

Ll=Li+2 1434 

M1=M1+1 1435 

M2=K2-1 1436 

GO TO 808 1437 

805 Jl=Jl+2 1438 

K1=K1+1 1439 

K2=K2+1 1440 

P1=M1+1 1441 
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M2=M2-1 1442 

GU TO 100 1443 

806 J2=J2+2 1444 

K1=K1-1 1445 

K2=K2+1 1446 

M1=MH-1 1447 

M2=K2tl 1448 

GU lU 100 1449 

807 Jl=Jl-2 1450 

K1=K1-1 1451 

K2=K2-1 1452 

M1=M1-1 1453 

M2=M^+1 1454 

GO ro 100 1455 

801 IX=JX+1 1456 

IF (IX-IS) bu9,8w9,8l0 1457 

809 IN0tX=l 1458 

GO TO 804 1459 

808 SuX=SCXt(2 .C=t=SPX) 1460 

GO TO 100 1461 

810 IV=iy+l 1462 

IF (lY-IS) 2811,2811,2812 1463 

2811 SUX=SJX-(2.C*( SPP-1.0)*SPX) 1464 

SOY=SQY+(2.C*SPy) 1465 

Jl=l 1466 

Kl=2 1467 

K2=K2+1 1468 

Li=2 1469 

L2=L2+2 1470 

Mi=l 1471 

K,2=M2 + 1 1472 

1X=1 1473 

811 II'J0hX=l 1474 

GO TO 100 1475 

2812 S0X=S0X+{2.C=i=SPX» 1476 

S0Y=S0Y-(2. C*( SPP-1 .0 ) *SPY ) 1477 

IF (JC-fMX + i) 698,851,851 1478 

851 K=1 1479 

SUX=SPX 1480 

S 0 Y= 6 t;Y+( 2 .C’i-'SPP=.'= 6 Py ) 1481 

1852 CJ.MINUE 1482 

1008 FU8F;AT (2110.5,15) 1483 

JC=1 1484 

GO 10 850 1485 

2226 IF (CC^JrC<N^^i)) 226 ,227,226 1486 

812 IF (LAST) 1001,1008,1001 1487 

1003 Cn.x'Tl.MUE 1488 

9941 CALL CALPLT < 24 .0, 0. 0,-3 ) 1489 

GU TG 9942 1490 

1001 WRITE (6,1002) I'^^l 

1002 FUKK.AT (5ri EiNjO ) 1^92 

CALL CALPLT (0.0,0.0,999) 1493 

RETURN l^*^^ 

END 
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Subroutine LINEl 


LINEl computes and stores all contour intercepts with side 1 in figure 2 for each 
pair of adjacent triangles in the 3 by 3 submatrix in turn. 

SUBROUTINE LINEl ( OEPTHL t DEP THM»C 1 , I NDEX , I tUl) 1496 

COMMON /3LK6/ SP 1497 

COMMON /BLK8/ SQX,SQY 1498 

COMMON /BLK9/ POS 1 X ( 2 0 ) , POS 1 Y ( 20 ) , CONT (20) 1499 

D1 =DE PTHL-UEPTHM 1500 

KX=OePTHL/Cl 1531 

SX=KX 1502 

IF (Dl.LT.O) CUNT (1 l=SX<'CH‘CI 1503 

IF (Di.GE.OJ CONF (1 )=SX*CI 1504 

104 CONTINUE 1505 

IF (Un l,2t2 1506 

1 CONTINUE 1507 

IF <CUNT( I)-D£PTHM.LE,0.0. ANQ.D1.lt. 0.0) 105,106 1508 

2 CONTINUE 1509 

IF ICCNT( I )-QEPrHM.LT.0.0.ANQ.Dl.GE.0.0) 106,105 1510 

105 IF (Dl.LT.O) DC=CONI (I )-DEPTHL 1511 

IF (Dl.GE.O) DC=QEPTHL-CGNT (1 ) 1512 

IF (Dt. £0.0.0) Ol=0.C0C000001 1513 

GO TO ( 1C7, 108, 109, 1 10 ), INDEX 1514 

107 POSlX(I)=SQX-( ( JC/AbS(Dl) )^SP) 1515 

POS1Y(I)-SOY 1516 

111 I-H-1 1517 

IF(I .GT. 20) PRINT 500 1518 

IF II.GT.20) B=(FUUR/0.0 )^>5^2 1519 

500 FORMAT (LHl^ CONTUUR IS lOU SMALL^/) 1520 

IF CDI.LT.0.0) CUNT (I )^C0NT (I -D+CI 1521 

IF (Ul.GE.0.0) CONT (1 )=CUNT (I -1)-CI 1522 

GO TO 104 1523 

1C8 POSlYd )=SUY+( (DC/ABS(D1) )^SP) 1524 

PUSlXiI)=SuX 1525 

GU TO 111 1526 

109 PJS1X( I ) = SUX+( (DC/AB$(D1) )^SP) 1527 

POSl Y( 1 )=S0Y 1528 

GU TU 111 1529 

110 PUSIX ( I ) = S0X 1530 

POSIY ( 1 ) = S0Y-( (UC/AbSIDl) ) «SP) 15 31 

GO TO 111 1532 

106 CGNTi 1 )=0.0 1533 

RETURN 1534 

END 1535 
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Subroutine LINE 2 


LINE2 computes and stores all contour intercepts with side 2 in figure 2 for each 
pair of adjacent triangles in the 3 by 3 submatrix in turn. 

SUBROUTINE LINE2 ( DE PTSM.UtPTHJ ,CI t I NDEX , 1 1 ,02 ) 1536 

CUPMCN /BLK6/ SP 1537 

COMMON /8LK8/ SOX.SOV 1538 

COMMON /BLK10/POS2X(20) ,P0S2Y(20) .CONTAl 20) 1539 

U2=DEPT HM-DEPTHJ 1540 

KX=DtPTHM/CI 1541 

SX=KX 1542 

IP (U 2 .LT. 0 ) CONrAIJJ )=SX*CI+CI 1543 

IP (O2.GE.0) CUNTAdl )=SX#Cl 1544 

115 CONTINUE 1545 

IP (02) 1,2,2 1546 

1 CONTINUE 1547 

IP ICUNTAI I I)-uEPTHJ.L£.O.AN0.02.LT.O) 116,117 1548 

2 CONTINUE 1549 

IF (CJNTAI I I )-DEPTMJ.LT.0.ANJ,U2.0£.C) 117,116. 1550 

116 IF (U2.LT.0) UC=CaNTA(II I-OEPTHM 1551 

IF (U2.GE.0) OC=DEPTHM-CUNIA(I I ) 1552 

IF (U2.Ea.0.0) 02 = 0.000000001 1553 

GO TO ( 1 18, 119 , 120, 121 ) , INDEX 1554 

llo POS2X( 1 I ) = S0X-SP 1555 

POS2YU I »=S0Y+( (UC/Ati6( 02) )*SP) 1556 

122 11=11+1 1557 

IF (11. GT. 20 ) 3 = (POUR/0.0)*=i'2 1558 

IF ( I I .GT .20) print 500 1559 

500 FORMAT!* CONTOUR INTERVAL IS TL<0 SMALL*/) 1560 

IF (U2.LT.O.O) CJNTAdI )=C0NTAdI -D+C1 1561 

IF (U2.GE.0.0) CUNTAdl )=C0.MTAdl -D-CI 1562 

GO TO 115 1563 

119 POS2X( d )=SUX+( (UC/A bS(02))*SP) 1564 

POS2Y( II )=SUY+SP 1565 

GO TO 122 1566 

120 P0S2X( d )=S0X+SP 1567 

POS2Y( d )=S0Y-( (0C/AGS(02))*SP) 1568 

GO TO 122 1569 

121 PoS2X( II )=S0X-( (0C/ABS(02) )*SP) 1570 

POS2Y ( d )=5GY-SP 1571 

GO TO 122 1572 

11 ? CONTA ( II )=0 .0 1573 

RETURN 1574 

END 1575 
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Subroutine LINES 


LINES computes and stores all contour intercepts with side S in figure 2 for each 
pair of adjacent triangles in the S by S submatrix in turn. 


SUBROUTINE LINE3 ( OEPTHL i JEP THJ ,CI . INOE X , 1 1 I 1 03 ) 1576 

COMMON /BLK6/ SP 1577 

COMMON /3LK.8/ SGX.SOY 1578 

COMMON /8LK11/PUS3X(20) ,PQS3Y(20» ,CONTB( 23» 1579 

03=DfcPTHL-DEPrHJ 1580 

KX=0EPTHL/C1 1581 

SX=KX 1582 

IF (D3.LT.0) CGNT8( I 1 n=SX«'C I+CI 1583 

IF ID3.GE.0) C0NT8( I I I )=SX*CI 1584 

341 CONTINUE 1585 

IF (03) 1,2,2 1586 

1 CONTINUE 1587 

IF (CONTBd I I )-OEPTHJ.LE,0 .0 .ANU.U3 .lt. 0.0)339, 340 1588 

2 CONTINUE 1589 

IF (CONTB( I II )-DEPThJ.LT.O.O.AND.03.GE.O .0) 340, 339 1590 

339 IF (03.LT.0) UC=CUNT B( i I I ) -OEPTHL 1591 

IF (U3.GE.0) 0C=UGPTmL-C0NT8( I I I) 1592 

IF (O3.E0.0.0) 03=C.OOC000001 1593 

GU TC (342,343,344,345) ,IN0EX 1594 

342 PUS3X( I II ) = S0X-( (UC/A8S( 03) ) -SP) 1595 

POS3Y( I 1 1 ) = SUY+( (OC/AB5( 03) ) ^i=SP ) 1596 

346 II1=IIH-1 1597 

IF{III.GT.20) PRINT 500 1598 

1F( 1 I I.GT.2G) 8 = ( FOUR/O.C ) <'<'2 1599 

500 FORMAT!* COCMiUUR INTERVAL IS TOO SMALL*/) 1600 

IF (U3.LT.0.0) CONTb( HI )=CGNT8( 1 I I-D+CI 1601 

IF (U3.GE.0.0) CJNT,3( I 1 I )=CUNT8( 11 l-D-Cl 1602 

GU TO 341 1603 

343 POS3XI I I I )=SUX+( (UC/ABS( 03) ) *SP) 1604 

PUS3Y( I I I ) = S0Y+ ( (UC/ABS( 03 ) ) *SP ) 1605 

GO TO 346 1606 

344 PUS3X( I I 1 ) = SUX+(0C/A8S( 03 ) )*SP 1607 

PUS3Y( HI )=S0Y-(UC/ABS(U3) )*SP 1608 

GU TO 34b 1609 

345 POS3X{ I I H=S0X-( ( 0C/ABS( C3) ) *SP) 1610 

PUS3Y( II H=SOY-( (OC/A8S( 03) ) *SP) 1611 

GU TO 346 1612 

340 CONTBI 1 1 1 ) = 0.0 1613 

RETURN 1614 

ENU 1615 
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Subroutine LINE4 

T.TN K4 computes and stores all contour intercepts with side 4 in figure 2 for each 
pair of adjacent triangles in the 3 by 3 submatrix in turn. 


SUbf<UUriNE Li;Mt4 (OEPTHJ.DfcPrHKtCl f lND£XtI4 ,D4) 1616 

CCPMON /BLK6/ SP 1617 

CUMMGN /8LKS/ SOX, SOY 1618 

COMMON /BLKi2/PUS4X(20),P0S4y(20),C0NrC(20) 1619 

U4=D£PTHJ-UZPTHK 1620 

KX=DtPrhJ/Ci 1621 

SX=KX 1622 

IF (D4.LT.0) COiMTC(l4 )=SX*CI+CI 1623 

IF (U4.Gt.O) CUNTCd'^ »=SX*Cl 1624 

207 CONTINUE 1625 

IF (04) 1.2,2 1626 

1 CONTINUE 1627 

IF (CONTC( I4)-i)tPTHK.LE. 0.0. ANU.D4 .lt. 0.0)208, 209 1628 

2 CONTINUE 1629 

IF (CONTC( I4)-0ePIHk. LI .C.O. AND. 04. GE. 0.0)209, .208 1630 

208 IF (D4.LT.0) UC=C0NTC(I4 )-UEPTHJ 1631 

IF (C4.Ub.O.O) DC = UEPTr)J-C3N rC { 14) 1632 

IF (04. EG. 0.0) U4=O.COCCCOC01 1633 

GU TO (210, 211 ,212,213) ,IN0dX 1634 

210 P0S4X(l4)=SUX-SP+( (UC/A0S(D4))*SP) 1635 

PUS4Y( I4) = SGY<-SP 1636 

214 14=14+1 1637 

1F(14 .GT.2C) PRINT 500 1638 

IF(I4.GT.20) 8=(FOUR/0.0)**2 1639 

500 FURMAT(« CONTOUR INTERV.AL IS TOO SMALL«'/) 1640 

IF (U4.LT.0.0) Cc)NTC( 14 )=C0NTC(I4 -D+Cl 1641 

IF (U4.GE.0.0) CuNTC(I4 )=C01TC(I4 -D-CI 1642 

GU TO 2C7 1643 

211 POS4X( l4)=SwX+SP 1644 

FUS4Y(14)=SGY+SP-( ( DC/ ABS ( 04 ) ) *SP) 1645 

GU TO 214 1646 

212 PUS4X( I4)=SQX+SP-( (UC/ABS(D4) )*SP) 1647 

PuS4Y ( 14 )=SGy-SP 1648 

GO TO 214 1649 

213 P0S4XI I4)=SWX-SP 1650 

POo4Y( l4)=S0V-SP+( ( 0 C/AdS( 04) ) »SP) 1651 

GO TO 214 1652 

209 CONTC( I4)=0.C 1653 

RETURN 1654 

END 1655 
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Subroutine LINES 


LINES computes and stores all contour intercepts with side S in figure 2 for each 
pair of adjacent trinities in the 3 by 3 submatrix in turn. 


SUaROUTINfc LIhj£5 ( DE PTriL f OEP FHK ,C 1 1 1 NDE X , I 5 ,0i>) 1656 

COMMON /8LK6/ SP 1657 

COMMON /BLK8/ SOX, SOY 1658 

COMMON /8LK13/POS5X(20) ,PQS5Y( 20J ,CONTO( 20) 1659 

U5=UEPTHL-DEPrHK 1660 

KX=UEPrHL/Cl 1661 

SX=KX 1662 

IF (U5.LI.0) CONIDU5 )=SX*CH-CI 1663 

IF (05.GE.0) CGNIUdb )=SX*CI 1664 

244 CONTINUE 1665 

IF (05) 1,2,2 1566 

1 CONTINUE 1667 

IF (C0NT0( I5)-DEPTHK.LE. 0.0. AND.05.lt. 0.0) 245,246 1668 

2 CONTINUE 1669 

IF (CONTU( I5)-OEPTHk.L1 .O.O.AND.D5.GE.O.O) 246,245 1570 

245 IF (05.LT.0) DC=C0NTD(I5 )-0EPTHL 1671 

IF (D5.GE.0.0) UC=OhPTHL-CJNTU ( 15) 1672 

IF (D5.E0.0.0) D5=O.OCOOC0001 1673 

GO TO (247,248,249,250) , INDEX 1674 

247 POS5XI 15 )=S0X 1675 

PUS5Y( 15)=S0Y<-( (0C/AbS(D5) )*SP) 1676 

251 15=15+1 1677 

IF(I5 .GT.20) PRINT 5CG 1678 

1F(15.GT.20) 8=(FUUk/0.0)*«'2 1679 

500 FORMAT!* CONTOUR INTERVAL IS TOO SMALL*/) 1680 

IF (D5.LT.0.0) CJNTDdS )=C0NTD(I5 -D+CI 1581 

IF (D5.GE.0.0) C0NTDd5 ) = CONTD(I5 -D-CI 1682 

GO TO 244 1683 

248 PUS5X(I5)=S0X+( (DC/AbS(D5) )*SP) 1684 

PUS5Y( I 5)=S0Y 1685 

GO TO 251 1686 

249 POS5X( I 5)=S0X 1687 

POS5Y(I5)=SOY-( ( JC/A3S(D5) )*SP) 1688 

GO TO 251 1689 

250 P0S5X( I5)=S0X-( (DC/AbS( 05) )*SP) 1690 

P065Y (I5)=SCY 1691 

GO TO 251 1692 

246 CONIC! 15)=0.0 1693 

RETURN 1694 

END 1695 
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Subroutine PLOT 


PLOT generates the necessary plotter instructions for drawing a straight line 
between two intercepts for a particular contour line. If the contour terminates or begins 
with a side which connects two boundary elements of the depth matrix, the contour value 
is plotted. 

PLOT is called from within a series of nested logical statements which first match 
the corresponding stored intercepts and then index sequentially until all contours within 
the confines of a triangle are plotted. 

SUHROUIINfc PLJT ( POS 3 X , POS 3 Y , P US 1 X , P OS 1 Y , PUS 2 X . P3S 2 V , CONT A , I S I 3N ) 1696 

CGMMUN /8LX7/ N . .'iN # P X 1 t P Y 1 , X MA X , YMA X t XMA Y f YMAY i Cl 1697 

CUMKUN /LIMITS/ M IN I ,MAX f , I AOOI 1698 

DATA CCT/U3770O00UCOCGC0C00000/ 1699 

CALL CAlPLF (PUS3X . P0S3Y ,3 ) 1700 

CALL CALPLT ( P JS 1 X , PUSl Y , 2 ) 1701 

DX = PtiS2Y-PY 1 1702 

DY=PCS2X-PX1 1703 

S=SU.<n (0X*DX) +DY^DY) 1704 

AMAX=XMAX+U .05 1705 

rfcMPX=CbNTA 1706 

IF (CCNTA.lt .103) GU TO 29 1707 

COM A=CuNTA- lAOUI 1708 

IF (S-U.07) 89,30,31.. 1709 

30 IF (Pno2X-').LJl5) 23,22,22 1710 

23 CALL COLUMN ( -C .0 d , P uS 2 Y , 0 . 0 7 , CuN T A , 0. 0 , -1 J 1711 

PY1=PCS2Y 1712 

PXi=PUS2X 1713 

GL T 0 29 1714 

22 IF (PUS2Y + 0.CJ 15) 24,24,25 1715 

25 CALL NUMOLk ( PC S2 X , C . C 5 , G . C 7 , CUNT A , 9 0. C , - 1 ) 1716 

PYi=POS2Y 1717 

PXi=P()S2X 1718 

GU TO 29 1719 

24IF(POS2Y + Y.MAX-J.0 515)26,27,27 1720 

2p call column ( PJS2X , YMAY , C. 07 ,CUNTA ,90.0 , -1 ) 1721 

PYi=PUS2Y 1722 

PXi=PQS2X 1723 

GU TO 29 1724 

27 IF (PUS2X-XMAX+0.0015) 29,29,28 1725 

28 call NUM3Ei\ ( AM AX , P0S2 Y , 0 . 0 7 , CUNTA , 0 .0 , - 1 ) 1726 

99 CGNfl.NUE 172 7 

PY1=P0S2Y 1728 

PX1=PUS2X 1729 

29 N=N-1*1SIGN 1730 

CnNTA=TEMPX 1731 

NN=KN+1 1732 

PE TURN 1733 

END 1734 
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Subroutine COLUMN 
By George Salley 


COLUMN right-justifies a BCD (binary-coded decimal) number. 

SUBRUUriME COLUMN ( X , Y , S t FPN tTH , N) 1735 

X AND y ARE THE CQORDlNArES OF THE L3WER RIGHTMOST EDGE 1736 

S IS THE SIZE OF CHARACTER TO 8E USED IN PRINTING 1737 

FRN IS AN ACTUAL FLOATING POINT NUMBER WHOSE YALUE IS TO BE 1738 

PRINTED ON THE PLOTTED OUTPUT 1739 

TH IS AN ANGLE (DEGREES! AT WHICil THE NUMBER IS TO APPEAR 1740 

N IS AN INTEGER SPECIFYING THE ACCURACY TO WHICH THE NUMBER IS 1741 

TO BE PRINTED 1742 

SPC = .857140 1743 

SPG =.285714 1744 

NL = 0 1745 

M = C 1746 

IF IN) 50,20,2U 1747 

20 NL = N+1 1748 

M = N 1749 

50 TFPN = RGUNC(FPN) 1750 

IFPN = TFPN*1J*«M 1751 

IF I IFPN ) 70,55 , /5 1752 

55 IF (IFPN) t)U,6u,100 1753 

60 NL = NL+1 1754 

GO TO 100 1755 

70 NL = NL+1 1756 

TFPN = AtiS ( TFPN ) 1757 

75 M = 0.4343*AL0G( TFPN)+1 .C 1758 

IF (M) 10U,10U,8j 1759 

80 NL = NL «• M 1760 

lOJ DLT =(SPC -S' S *1 FLJAl ( NL ) ) )- ( SPG*S ) 1761 

I = TH*0.C17453 1762 

YP = Y -( DL I *S IN ( T ) ) 1 763 

XP = X -(OLT'-'CUBm ) 1764 

120 CALL NUMBER ( XP , YP , S , FPN , I H , N ) 1765 

RETURN 1760 

END 1767 
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OPERATING INSTRUCTIONS 


Input 

Basically the input data can be divided into two categories. One category is 
termed "program data," and the other is referred to as "control-point data." (The max- 
imum number of control points is limited to 1000.) Program data are composed of first 
a format (xiser's choice) for reading the floatir^-point program variables; the fixed-point 
program variables and the appropriate codes for each option follow. If the control-point 
data are not on tape, the user must supply a format for reading each control-point data 
card. The control-point data are composed of individual cards or tape records each of 
which contains the X-Y coordinate of each scalar variation and its magnitude Z, in 
that order. 

The control-point data can originate from any source. If an algebraic expression 
is to be utilized to obtain the control-point data, irregularly spaced scalar variations are 
preferred. Hov/ever, regularly spaced scalar variations can be processed with the 
program. 

In view of the various options available to the user (to be discussed subsequently), 
the number of required data cards is not rigidly fixed. Thus, a figure indicating the card 
on which a particular variable is to be found is not possible. However, in figure 4 the 
order in which each variable is read and the format utilized for reading it are indicated. 

In table n program variables are listed and its function in the program is defined. 

Some knowlecfee of the input data is required to adequately assign the program 
variable values. Caution should be exercised in the selection of DIST (grid separation), 
since it is primarily responsible for establishing the dimensions of the depth matrix. If 
DIST is too small, the maximum dimensions of the depth matrix (60 by 60) will be 
exceeded. DIST also controls the dimensions of the region within which control points 
can be found. (See the discussion on "Depth Matrix.") All the control-point data are 
searched for values which fall within a ring with the initial inner radius set equal to zero 
and the outer radius set equal to the incremental radial difference of DIST/lO. The ring 
has a mean radius of one -half the sum of the two radii. Since the mean radius of the ring 
is incremented in multiples of the incremental radial difference before each search of the 
data, DIST probably influences the required computer time more than any of the remaining 
variables. 

The imits for DIST, X, Y, and DISPGRD are required to be the same. If X 
and Y represent different quantities, such as Mach number and height, respectively, 
an appropriate scale factor must be applied to X to establish a fictitious one-to-one 
correspondence between the magnitudes of the two variables. Since the objective is to 
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establish roughly the same order of magnitude between the two variables, Mach number 
should be converted into the units of height per second. 

It is assumed that X, Y, and sometimes Z are expressed as integral multiples 
of the user’s basic unit (meters, inches, etc.). Thus, if it is desired to contour data 
which are expressible in fractional parts of the basic unit, X, Y, and Z must be 
scaled up in magnitude so that the significant part thereof is plotted. That is, data 
expressed to 0.0001 inch are to be multiplied by 10 000. If similarly treated, SCLE will 
reestablish the desired plot dimensions. 

The inputs MAXI, MINI, TANG, and IDMX establish the criteria for either accepting 
or rejecting a given control point. The upper and lower limits for the magnitudes of the 
control -point X-component are set by MAXI and MINI. If both MAXI and MINI are zero, 
no data points are rejected. Together TANG and IDMX determine whether the relative 
height difference between the Z -coordinates of two control points is sufficiently small to 
permit their combined use in the gradient determinations. Of the two variables TANG 
is the dominant variable over a radial distauce equal to the product of TANG and IDMX. 
Thereafter, IDMX becomes the dominant variable. The required computer time can be 
affected by TANG and IDMX. 

The units for Z, IDMX, and Cl are the same. The input Z may be any scalar 
variation (height, temperature, etc.). The program assumes IDMX, Cl, and Z are 
integral multiples of a basic imit and not fractional parts thereof. 

Options 

Control-point data . - The control-point data may be on either punched cards or tape 
and may be either floating point or fixed point. Program codes specify which of the 
four options are to be expected. If the punched-card option is selected, a format for 
reading each data card must be furnished by the user. The subroutine READIT will 
automatically place the data in the proper form and establish the proper t}q)e required 
for processing. 

Sorting of control-point data .- There is no requirement for the control-point data 
to be in a specific order. The program subroutine GRIDIT will automatically scan all 
the data furnished until a maximum of KMAX control points are selected for the gradient 
computation. The scanning process is, however, a very time-consuming activity. In 
fact, it is by far the largest user of computer time and should be minimized as much as 
possible. As a result, the program contains an optional sortii^ routine, SORTIT. 

SORTIT can be called with a code of +1 for ISRT. The variable ISRTOPT specifies the 
manner in which the data are to be sorted. 
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Interpolati on.- The digital program will expand each 3 by 3 submatrix into a 
2SPP + 1 by 2SPP + 1 submatrix in AUTOCON. Each missing element, if any, is deter- 
mined by use of linear interpolation. The input SPP permits the size of each triangle 
over which the line segments, representing contour lines, are drawn to be reduced so 
that the contours will appear more continuous. The value of SPP can not exceed 8. 

Deck Setup 

As stated, the digital contouring program is composed of seven overlaid programs 
and their various subroutines. CONTOUR is the main overlay (0,0 level) and is, as a 
result, the program identification. 

A deck setup for a case in which the control-point data are contained on a deck of 
pimched cards is shown in figure 5, and a deck setup for a case in which a binary tape file 
is used is shown in figure 6. The binary tape must be one which was generated by using 
RECOUT. If the control-point data are contained on a deck of punched cards, the control 
cards for tape file 9 are not needed. 


Output 

The digital program outputs data in two forms. First, there is a printed listing of 
information pertaining to the depth matrix (dimensions, origin, etc.) and to the plot size, 
as well as a frequency distribution plot of all the differences between the actual scalar 
measurement at each control point and that predicted by the depth matrix. Second, 
plotter instructions for a graphic display of the scalar variations are written on tape. 
The tape is then used by the mechanical plotter to produce the contour charts. 

Several diagnostics have been programed into the digital contouring program. 

Each attempts to isolate potential execution difficulties and to inform the user before 
they occur, in the interest of saving computer time. The following reasons are given 
for the diagnostics: 

Control-point data exceed 1000 points. 

Maximum dimensions of depth matrix are exceeded. 

Depth-matrix dimensions are too small. 

Minimum depth- matrix elevation is too small. 

Depth matrix cannot be smoothed. 

Contour interval is too small. 

A mode 4 stop is anticipated for each diagnostic. The diagnostics themselves are some- 
times lengthy and are not given herein. 
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JOB, 1,1000,^000. — 

RUN(S) 

NORFL. 

SETIHDF . 

LGO. 

REWIND (CALTPE) 

REQUEST TAPEXX,HI,X. 
REWIND (TAPE XX) 

COPYBF ( CALTPE , TAPE XX) 
UNLOAD (TAPE XX) 

EXIT. 

UNLOAD (TAPE XX) 

7,8,9 

Source Deck 

7.8.9 

Program Data 
Control Point Data 

6.7. 8.9 

NOTE: TAPEXX may be 

Figure 5-- Deck setup 
of control-point 


JOB, 1,1000, ^000. 

RUW(S) 

NORFL. 

SETINDF . 

REQUEST TAPE1,HY . XXXXXX ,R0L, 
REWIND (TAPEl) 

CALTP , RIL, REWIND (TAPE9) 

COPYBF (TAPEl, TAPE9) 

REWIND (TAPE9) 

DROPFIL(TAPEl) 

LGO. 

REWlND( CALTPE) 

REQUEST TAPEXX, HI,X. CALTP , RIL, 

REWIND( TAPEXX) 

COPYBF ( CALTPE, TAPEXX) 

UNLOAD(TAPEXX) 

EXIT 

UNLQAD(TAPEXX) 

7,8,9 


Source Deck 


any legal file name. 

when a punched deck 
data is used. 

7 . 8.9 

Program Data 

6 . 7 . 8. 9 


NOTE: TAPEXX may be any legal file name 

XXXXXX may be any legal tape number 

Figiire 6.- Deck setup when a binaiy file of 
control-point data is used. 
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II III 


SAMPLE PROBLEM 

A sample problem was set up for a punched-card deck of control-point data. For 
each program data card, the magnitude of each variable punched on each data card is 
shown in figure 7, and in table HI, the control-point data read by the variable format 
FRMT ((3110) in fig. 7) are tabulated. The plotted graphic is shown in figure 8. 


Card 

Column number 

type 

123^567890123^567890123^5678901234567890123^^567890123^56789012345678901234567890 

1 

( 4 ei 8 . 8 / 4 ei 8 . 8 ) 

2 

10 10 6000 

5 

+ 30000000 E +04 

+ 12000000 E +02 + 10000000 E +01 

+ 10000000 E +01 

k 

+IOOOOOOOE+O3 

+9OOOOOOOE+O5 +I2OOOOOOE+O5 

+ 10000000 E +01 

5 

-1 -3600 

9000 


6 

-1 +1 +0 



7 

(3110) 




Figure 7 .- Program data for sample problem. 


The computer program requires approximately 65 000 octal words of storage, and 
the execution time varies from job to job. The execution time is primarily dependent on 
the number of control points being processed, the distribution of the control points, and 
the magnitude of DIST. On the average, 1-| seconds of central processor time and an 
eqvlivalent amount of peripheral time is required for each control point being processed. 

It is possible to reduce the time required to process a particular job by increasing DIST. 
The sample problem required approximately 739.95 decimal seconds of central processor 
time and 767.18 decimal seconds of peripheral time. The peripheral time may vary con- 
siderably even if the same data are used. 
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Listing of Unsmoothed Depth-Matrix Data 


THE FOLLOWIN.5 IS INFORMATION PERTAINING TO THE DEPTH MATR I X— S IZ E , D I MENS I ONS , ETC . 


CORNER COORDINATES OF AREA TO BE CONTOURED 

NORTH WEST CORNER — X= A. OCOOOOOOE + 03 Y= 6, 00307200E 
SOUTH EAST CORNER — X= 5 . 098 80000E+04 Y= 0. 


DIMENSIONS OF REQUIRED PLOTTING SURFACE EXCLUSIVE OF GRID 

Y DIMENSION OF DEPTH MATRIX IS 8 . 00400 I NCHE S 

X DIMENSION OF DEPTH MATRIX IS 6.26493 INCHES 

22ELEMENTS IN Y DIRECTION OF DEPTH MATRIX 
17ELEMENTS IN X DIRECTION OF DEPTH MATRIX 

YOU HAVE ELECTED TO REJECT ALL CONTROL POINTS WITH SCALAR VARIATIONS LESS THAN -3600 AND GREATER THAN 
THUS, 

YOU HAVE A MINIMUM SCALAR VARIATION OF -3554 

YOU HAVE A MAXIMUM SCALAR VARIATION OF 327 


THE FOLLOWING INFORMATION REGARDING YOUR UNSMOOTHED DEPTH MATRIX IS FURNISHED 


ADJUSTED BOUNDARIES OF DEPTH MATRIX 

MINIMUM Y 0. MAXIMUM Y 6. 003 07200E+04 

MINIMUM X 4.00000000E+03 MAXIMUM X 5 , 09880000E+04 


UNSMOOTHED DEPTH MATRIX 


946 1 

9000 

8504 

7903 

0021 

8123 

7962 

8055 

8 107 

7914 

7794 

7637 

7612 

7586 

7592 

7784 

7953 

9457 

9010 

8559 

794 4 

8029 

8185 

8082 

8188 

8247 

8105 

8042 

7886 

7947 

7960 

7979 

8232 

8130 

9463 

901 5 

8575 

8009 

R098 

8257 

8209 

8339 

8383 

8333 

8299 

7873 

8622 

8201 

8 409 

8391 

8311 

9452 

9030 

8601 

8137 

8085 

8313 

8363 

8422 

8642 

8551 

8559 

8403 

8753 

8538 

8606 

85 81 

8529 

9440 

9045 

9601 

B25 0 

8210 

84 56 

0529 

8643 

8797 

8829 

8745 

8677 

8876 

8934 

8703 

8695 

8628 

9458 

9052 

8656 

8370 

8135 

8639 

8726 

8745 

9004 

9155 

8956 

905 1 

8970 

8929 

8900 

8774 

8681 

947 0 

9060 

8687 

834 1 

8442 

86 82 

8934 

9116 

9117 

9263 

9177 

8925 

9413 

9042 

8956 

8888 

0 793 

9460 

9093 

8729 

8468 

8 546 

8898 

9177 

9326 

9 340 

9306 

9326 

9232 

9277 

9242 

9074 

8902 

8806 

9492 

9077 

8785 

860 8 

8465 

92 35 

9404 

9499 

9411 

9542 

9448 

9490 

9313 

9231 

9174 

8950 

86 70 

9461 

9131 

8802 

8682 

8969 

9149 

9531 

9562 

9631 

9525 

9563 

9504 

9237 

9332 

9113 

8716 

8772 

9482 

9114 

8836 

8786 

9178 

9282 

9 775 

9668 

9672 

9540 

9568 

9423 

9487 

9070 

8808 

8482 

8489 

9450 

9163 

8839 

6918 

9332 

9662 

9 756 

9731 

9770 

9512 

9507 

9437 

9172 

8891 

8556 

8171 

7841 

9453 

9248 

8 804 

9052 

9403 

9747 

10 026 

9778 

9596 

9672 

9532 

9284 

9013 

8781 

8102 

7821 

7196 

9489 

9122 

9016 

9180 

9859 

9840 

9942 

9813 

9699 

9243 

9823 

9159 

8916 

8323 

7907 

7357 

6767 

9478 

9297 

8936 

9341 

9881 

9985 

9996 

9805 

9720 

9566 

9291 

8968 

8568 

8076 

7520 

6710 

6215 

9499 

9245 

9076 

9539 

9969 

10010 

10072 

9832 

9639 

9573 

9067 

8730 

8237 

7701 

6909 

6088 

5572 

9485 

9258 

9178 

9 791 

10048 

10136 

10024 

9827 

9674 

9283 

9015 

8582 

7800 

7163 

6339 

5534 

4903 

9529 

9371 

92 90 

9879 

10115 

10137 

10070 

9817 

9491 

9261 

8776 

8341 

7507 

6738 

5753 

4953 

4217 

9517 

9277 

9512 

10000 

10221 

10163 

10125 

9751 

9483 

9046 

8651 

7912 

7033 

6186 

5239 

4379 

3568 

9517 

9460 

9449 

1018 1 

10303 

10192 

10111 

9727 

9345 

8886 

8537 

7493 

6576 

5624 

4726 

3814 

2958 

9556 

9398 

9782 

10107 

10322 

10231 

101 13 

96 58 

9263 

8769 

8214 

7271 

6071 

5068 

4139 

3224 

2338 

950 7 

9351 

9790 

10188 

10318 

10180 

10057 

9609 

9195 

0609 

7911 

7132 

5660 

4571 

3605 

2662 

1736 


9000 



Listing of Smoothed Depth-Matrix Data 


THE FOLinwINR INFORMiTICN PcnARDINR T^E S^mTHEH DEPTH M^TRIX IS FUR^1ISHE^ 


IT HAS RPFN nETERMI'^iFr) T«AT CACH ELEMENT TF TOUR UNSMnOTHED DEPTH MATRIX IS INAnEQUATELY FILLED 


ADJUSTED RD1JNUART=S 1= D'‘PTH MATRIX 

MINIMUM Y T. MAXIMUM Y 6 .00 S'*? 2^ ''E +T4 

MINIMUM X JC + -3 maximum x 5 • C 9 P8f C 00F + D4 


REf;PTD WILL NOT PFRMTT TMF INTroNAL ELEMENTS OF I TO EXCEED THE RDLLOWINO LIMITS 

MINIMUM FIFVATION PRIIR TO SMOOTHINC = -2^2/* 

MAXIMUM FLFVATION PRIDR TH SMOPTHTNr, = 322 

RE PAUTiniS OF ELFMFNT 

MINIMUM ELEVATION AFTER I NTFB POL AT ION FTC = -3424 

MAXIMUM FXPECT«=0 SLOPF RFTWFFN TWO CONSECUTIVE r7RID POINTS = l.ODTCO 


SMOOTHED DEPTH MATRIX 


9461 

OOOO 

P534 

7933 

83 21 

R123 

7962 

8355 

8107 

7914 

7794 

7637 

7612 

7586 

7592 

7784 

7953 

94'j-' 

orin 

8569 

39 44 

8T?c 

PIPK 

pro 7 

8188 

R247 

P10 5 

8042 

7806 

7947 

7560 

7979 

8232 

8130 

946 3 

9015 

R575 

83-^0 

83QR 

P7K7 

8209 

0339 

8383 

8333 

8299 

7873 

8622 

9201 

843 9 

P391 

8 311 

94*^2 

OP^A 

R6m 

8137 

8385 

P31 3 

R36 3 

8422 

8642 

8551 

8559 

8403 

8753 

8538 

8606 

P561 

8529 

94^3 

9P45 

8601 

8253 

871 r 

8456 

0529 

P643 

R797 

R829 

874 = 

8677 

8876 

8534 

8 70? 

8695 

862P 

94SR 

9P52 

86 58 

P37"> 

Q 1 85 

P630 

8726 

8745 

9004 

0155 

8956 

9051 

R97C 

8929 

890C 

8774 

8681 

947'> 

9 063 

8687 

8341 

8442 

P68? 

R93A 

9116 

9117 

0263 

9177 

8925 

9413 

9C42 

8956 

8888 

8793 

946'' 

9093 

8739 

8468 

QS4f 

RRoa 

9177 

9326 

934'* 

93 16 

9326 

9232 

9277 

9242 

9074 

8902 

P806 

9403 

9077 

87R5 

86)8 

846^ 

9238 

9404 

94cq 

9411 

9542 

044 8 

949') 

9313 

9231 

9174 

80 = 0 

0673 

9A.t I 

9131 

RRP3 

8 68 7 

8960 

9149 

9 = 31 

056? 

96 31 

9=25 

9563 

0534 

9237 

9332 

9113 

8716 

8772 

9^‘^RP 

0114 

8836 

8 7 8.6 

Q17P 

Q7R 2 

0775 

9668 

9672 

9 540 

956R 

9423 

94P7 

9C70 

8 9C8 

P402 

8409 

94R3 

9163 

8830 

8918 

933? 

Of 62 

9756 

97 ? 1 

9773 

951 2 

0637 

9437 

9172 

8 891 

8556 

8171 

7841 

94F3 

934B 

RR6 4 

0052 

04''3 

9747 

10026 

0778 

9596 

9672 

9532 

9284 

9013 

8781 

8132 

7821 

7196 

94BO 

913? 

9016 

9180 

9 8 59 

9840 

9947 

9 813 

9609 

9243 

9823 

9159 

8916 

8323 

7987 

7357 

6767 

9470 

9 297 

R93R 

0 3 VI 

OR ni 

9988 

9096 

9 53 8 

9723 

0566 

0291 

8968 

P56R 

8C76 

7523 

6710 

6738 

RA-Oo 

o?45 

90 7 R 

9539 

9069 

10310 

10072 

03?2 

9639 

9573 

9067 

8730 

8237 

7701 

6939 

6P09 

6709 

940*; 

925R 

91 78 

9791 

IP '>48 

13136 

nr 24 

9827 

9674 

92R3 

901 5 

8582 

7883 

7163 

7036 

6922 

68 OR 

9S29 

9 371 

9790 

9870 

101 15 

' 3137 

1C 370 

C8l- 

9491 

9261 

0776 

8341 

7507 

6738 

6887 

69C4 

6«21 

9517 

9 277 

9512 

1 9333 

V'221 

1D163 

10125 

97=1 

9483 

9346 

8651 

7912 

7033 

6P85 

&RR6 

6805 

6904 

9517 

046'! 

9449 

131 81 

13303 

iri92 

K 111 

0727 

9345 

88B6 

8537 

7493 

6576 

6730 

6808 

6fl51 

6 894 

955 6 

939R 

97B2 

13137 

13322 

10731 

K113 

965 n 

9263 

876 9 

8214 

7271 

6923 

6826 

6817 

6 8 34 

6851 

95D7 

9351 

9793 

13188 

13318 

131R0 

10C57 

9 63 9 

9195 

R6f 9 

7011 

7132 

7270 

6522 

6826 

6 8L.7 

6008 


YOUR SMTOTHEO DEPTH MATRIX. I. IS PRINTED ABOVE. 
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Sample Problem Output for Evaluation of Smoothed Depth Matrix 


*T T4IS PnrNT THF SC AL VfRUTlONS PREOICTED BY I THE LOCATION OF FACH 
IS rnMPARFn with the SCAI&P MACNITUnE OF THF CCNTROL POINT. 

^ PLOT FOR FREOUFNCV OF OrCURRENCF VS DISCREPANCY FOLLOWS 4S PART OF THIS L 


CONTRCL POINT 
I STING 


FREQUENCY OF OCCURRENCE 

■>1 734S67 

n 

I -SQt* 

f;; 

p 7 + * + + + * + 

P ♦*«*](!«**%*♦♦♦* ♦;<«♦*«♦« 

a + + + 

rg + 

r 9^* 

Y 119 * 

1 39 * 

ALL CONTROJ POINTS HAVf RF*=N CONPAPFQ AND THE RESULTS AR£ SHOWN 
IT WAS n^TERMINFO FRQvi T COMPARISONS THAT T1<FPE WAS A 
MFAN QTFFERENCF = 7 . ' S'* c-r 

AND A 

ST OrV = ? . 1 1 9^1 S 7 ?F- -iT 

THF OFOTH MATRIX WILL NOW PLOTTEO 

A TWO niWFNSTONAL CONTOUR CHART WILL BF PRFRAREP WITH A SUITABLF GRID 

EACH GRin LINF represents AN iNT^nRAL MULTIPLE OF OISPORQISEF INPUT INSTRUCTIONS) 

AND THF CONTOURS API RnsiTTONFO RFLATIVF TO THIS GRID 

PIOTTTNG RFGINS 'WITi T *F IJOPER LEFT OANO CCRN'FR OF THE DfPTH MATRIX 


P 


9 lOX 30 number 

5 

6 
24 
51 

120 

223 

88 

33 

16 

4 

5 
4 


FNO 
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FORTRAN VARIABLE NAMES 


The major FORTRAN variable names and arrays are defined in table IV. Where 
required, the algebraic equivalent of the variable is shown. Certain variables assumed 
different names in different programs and subroutines and are so indicated in the table. 

CONCLUDING REMARKS 

A generalized digital contouring program has been developed and described. The 
program will accept any type of data regardless of the units or the dimensions of the 
region over which the data are to be plotted. Further, the program makes no assump- 
tions regarding either the absolute or the relative magnitudes of the scalar variations 
being plotted. 

The digital contouring program is limited in the number of acceptable control 
points and, to a degree, in the magnitude of the grid separation. A maximum of 
1000 control points can be processed at one time. The grid separation should not be so 
small that it requires more than 60 grid points in any one direction in the depth matrix 
which represents the region being contoured. Further, the depth matrix cannot be 
smaller than a 3 by 3 array. 

A sample problem with 552 data points was processed. The resulting computer 
listing of the output is shown, as well as the machine -generated contour chart. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., September 3, 1970. 
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TABLE I.- FUNCTION AND STORAGE OF EACH OVERLAID PROGRAM AND SUBROUTINE 


Overlay 

Program 

Subroutine 

Storage 

Description 

Page 

(0, 0) 



16 460q 


13 


CONTOUR 


12 327g 

Reads program options and directs the order 
of execution of the overlays 


(1,0) 



5 625g 


14 


READIT 


442g 

Reads either fixed-point or floating-point 
X-, Y-, and Z- coordinates from either 
card or tape and generates a coded 
binary tape of data suitable for use by 
the remaining overlays and subroutines 




RECm 

302g 

A system routine for reading binary and 
BCD tapes generated by RECOUT 




RECOUT 

340g 

A system routine for packing binary and/or 
BCD shorter records into longer logical 
records of approximately 512 jq core 
memory words 


(2, 0) 



27 305g 


16 


SORTIT 


164g 

Prepares the input to SORT2 




SORT2 

207g 

A routine for calling SORT/MERGE using 
RECIN and RECOUT (Data is ordered 
according to program options.) 


(3, 0) 



13 060g 


17 


COUPLE 


13 060g 

Manages lower level overlays 


(3, 1) 



31 712g 


17 


MATCAL 


23g 

Manages the subroutines GRIDIT 
and REGRID 




GRIDIT 

17 160g 

Evaluates the gradient at each control 
point and evaluates the depth matrix 

18 



REGRID 

2 260g 

Checks each element of the depth matrix; 
fills in the missing elements; smooths 
the matrix 

25 



CONVERT 

168 

Converts the units of the input data to the 
plot equivalent in inches 

33 



TESTI 

251g 

Examines the boundary elements of the 
depth matrix and reduces the matrix 
dimensions when all the elements in a 
column or row are missing 

31 



CKPOINT 

107g 

Checks data for control points which have 
been observed more than once 

30 

(3, 3) 



27 061g 


34 


CKZFIT 


12 634g 

Compares each scalar variation predicted 
by the depth matrix at the location of 
each control point with its control-point 
counterpart 




DISCOT 

555g 

Interpolates between depth-matrix elements 
to evaluate the predicted scalar variation 
at the location ^ each control point 




SUBl 

5 155g 

Compares the predicted scalar variation 
with its control-point counterpart and 
plots a frequency distribution curve on 
the listing 

36 
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TABLE I.- FUNCTION AND STORAGE OF EACH OVERLAID PROGRAM AND SUBROUTINE - Concluded 


Overlay 

Program 

Subroutine 

Storage 

Description 

Page 

(3, 2) 



25 056g 


39 


AUTOCON 


3 351g 

Determines grid spacing and plots origin, as well 
as controls the orderly processing of the columns 
and rows of depth matrix and controls the orderly 
call to each subroutine 




GRID 

117g 

A CalComp modular routine for plotting a grid of 
equally spaced lines 




CALPLT 

154g 

Basic CalComp instruction to plotter pen; estab- 
lishes plot origin and pen position (up or down) 
before, during, and after changing locations 




LINEl 

223g 

Determines and stores the coordinates of the inter- 
cept of each contour line with side 1 of a triangle 

47 



LINE 2 

226g 

Determines and stores the coordinates of the inter- 
cept of each contour line with side 2 of a triangle 

48 



LINE 3 

245g 

Determines and stores the coordinates of the inter- 
cept of each contour line with the common side 
(side 3) of adjacent triangles 

49 



LINE4 

230g 

Determines and stores the coordinates of the inter- 
cept of each contour line with side 4 of the 
adjacent triangle 

50 



LINES 

223g 

Determines and stores the coordinates of the inter- 
cept of each contour line with side 5 of the 
adjacent triangle 

51 



PLOT 

1 366g 

Plots each line segment connecting corresponding 
contour intercepts on opposite sides of a triangle 
and calls NUMBER 

52 



NUMBER 

246g 

Labels those contour lines which terminate at the 
boundaries of the plot 




LEROY 

25g 

Slows plotter for Leroy pen 




CALCOMP 

1 056g 

Selects CalComp plotting system 




COLUMN 

163g 

Right-justifies BCD number 

53 



AXES 

510g 

Plots and annotates axes 
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TABLE II.- LIST OF VARIABLE NAMES AND DEFINITIONS 


FORTRAN 

variable 

name 

Card 

type 

Definition 

FRMT 

1 

Variable format for reading data on card 
types 3 and 4 

KMAX 

2 

Integer specifying the number of control- 
point gradients required to compute an 
elevation at a particular grid point 

KMAXl 

2 

Integer specifying the number of control 
points required to compute a gradient at 
a particular control point 

IDMX 

2 

Integer specifying the maximum difference 
between scalar variations 

DIST 

3 

Floating-point variable specifying the spacing 
between elements of the depth matrix 

UNIT 

3 

Floating-point variable for converting the 
units of DIST into the smallest unit for 
which a conversion factor is available 
for converting UNIT to inches 

CNVTOIN 

3 

Floating-point conversion factor for con- 
verting UNIT to inches 

SPP 

3 

Floating-point interpolation factor for sub- 
dividing the submatrices in AUTOCON 

Cl 

4 

Contour interval 

SOLE 

4 

Scale of the graphic 

DISPGRD 

4 

Desired separation between grid lines of the 
graphic 

TANG 

4 

Maximum anticipated slope between control 
points 

ISRT 

5 

Code specifying whether data are to be sorted 



+1 sort data 

-1 do not sort data 

MINI 

5 

Minimum anticipated scalar variation 

MAXI 

5 

Maximum anticipated scalar variation 

IRDTP 

6 

A code specifying whether control-point data 
are on tape or punched cards 



+1 data are on tape file 9 

-1 data are on punched cards 
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TABLE II.- LIST OF VARIABLE NAMES AND DEFINITIONS - Concluded 


FORTRAN 

variable 

name 


Card 

type 


IPTOPT 


6 


ISRTOPT 


6 


FRMT 


7 


Definition 


A code indicating type of data 

+1 data are fixed point 

-1 data are floating point 

A code indicating primary and secondary 
variables in SORT/MERGE routines 

Code Variable Direction Order 

-1 X Ascending Secondary 

Y Descending Primary 

0 X Ascending Secondary 

IDIST Descending Primary 

1 IDIST Ascending Primary 

Y Descending Secondary 

2 Y Descending Secondary 

X Ascending Primary 

A variable format for reading X, Y, 

and Z parameters on control-point data 
cards 
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TABLE m.- TABULATED CONTROL POINT DATA 
FOR SAMPLE PROBLEM 


IX 

lY 


CC 

lY 

JZ 

^onr 

* ’0 

-445 

676 3 

564 2 3 


/.inn 

1 2 '"‘C 

-460 

6 90 3 

57629 

-9P.1 


2401 

-471 

691 9 

5R829 

-9P6 


36'“'2 

-479 

693 9 

600 7 9 

-902 


4RA? 

-493 

7979 

n 

-*=^08 

43?7 

6302 

-504 

7987 

12’''0 


4^4*7 

7203 

-511 

7999 

24'i? 

-571 

^387 

B4''4 

-520 

3079 

36 0 3 

-S97 

439^ 

960 S 

-576 

373 c 

4803 

-61 9 

440 7 


-53? 

3 33 9 

60 0 3 

— 644 

44 R 7 

1279? 


3399 

720^:. 

-67^ 

5.f47 

13 206 

-565 

341 1 

3405 

-69P 

4747 

14407 

-^30 

3551 

960 5 

-*^22 

4R0 7 

15607 

— *;01 

8*^67 

1 

-75P 

/ qi c 

16337 

-596 

P537 

1 2 6 

-7-^7 

A7''> q 

iRfV^Pt 

-61 6 

3627 

132' 7 

-303 

407 c 

19 709 

-620 

363'^ 

14 408 

-3?P 

400K 

2041 n 

-6 26 

3647 

156'*' 3 

-331 


2161 

-63*“ 

3777 

168 ')9 

-R7S 

^ -4 3 

22R11 


8337 

iRino 

-Q 


24.712 

-646 

39P7 

1921 0 

-09R 

*51 3 « 

2571 3 

-6-9 

9047 

2 0 <• 1 ' 

— r 

^795 

2641 3 

-630 

qOSC 

2161 0 

— OTf 

5705 

276^ 3 

- 6 OR 

0 1C 9 

228 1 I 

-1 on/- 

iS4*5 c 

23314 

-701 

^21 5 

2401 3 

-1 03? 

54^7 

30015 

-7*^3 

9735 

252** 3 

-^*1 oq/.. 

56<^7 

3121 6 

-7 7? 

977S 

2641 3 

- 1 ^ 6 0 

^S73 

3 2416 


Q? P ? 

2 7M4 

_ 1 P 3 

*^64 3 

33616 

-729 

0 ?0 R 

2 3 3 


^683 

34 317 

-73’"' 

CRTS 

3 "16 

- 1 1 0 

S693 

3601 R 


0 63 0 

31 21 6 

- ’ 13^ 

*5 70 3 

37 7 1 R 

-74 if, 

963’*^ 


- 1 *1 r T 

^7B3 

3341 3 

-760 

969’^ 

3361 7 

-1 17*^ 

so/. 3 

39619 

-737 

070 7 

3-^ 81 8 

- 1 1 3 3 

604 3 

403 20 

-- p n r 

934-7 

36"' 10 

- 21^ 

6103 

47 021 

-319 

936 3 

3771 9 

2^2 

411c 

43771 

- 3 2 5 

9333 

38419 

-1 257 

62*5 5 

444 7 2 

-BCD 

QO?3 

3967 0 

-1274 

6 77 1 

4S673 

-3^'6 

QQR1 

4f 321 

-I29r 

6^91 

46823 

-363 

9943 

47 071 

-*» 3^7 

6331 

480 2 3 

-37? 

13023 

43221 

-1331 

633 S 

492 2 3 

-3 76 

1 0 1 a 3 

4 4 7 7 

-1 36S 

6 3? 1 

504? 

-*^31 

1,0 23 2 

45 S? 3 

-1390 

643 J 

SI 6 26 

-396 

V*' 7 3 

46 8 74 

-141 1 

6691 

52326 

-02? 

1 03*^*^ 

48076 

-1 474 

A691 

S40 27 

— 04 r. 

1 3 49 5 

4922s 

-16PS 

67?1 

55227 

-962 

1.C51 1 

5C 476 

-1 A 69 




*» 0 5 3 1 

51626 

-1 68? 


TABLE m.- TABULATED CONTROL POINT DATA 
FOR SAMPLE PROBLEM - Continued 


DC 



DC 

rr 

IZ 


5732? 

-1490 

1457 3 

50476 

-1 871 

1 r 7 9 

540 2? 

51 0 

14 5 8 2 

51676 

-1903 

1 ^. c > \ 

55727 

-1527 

1 4*=^95 

67826 

-1^30 

1 *^^67 1 

5^,479 

-1 544 

14735 

647 77 

-1<767 

1 

57f 37 

57g- 

147*= \ 

5527R 

-1093 

1 r q*5 1 

S R R ^ '*:- 

-I603 

14771 

66478 

-20! 8 

1 

«S43 3”) 

-1671 

1481 1 

57628 

- ?n/->3 

1 3 


<5? 

1 4 8 1 9 

6887P 

-2 06 8 

^214? 

1 2 “ 1 

77 

14831 

6rr 29 

-2^90 


2 4- 1 

/r 

l*^^l 1 


377 

T 2 1 T q 

335 1 

7 

! 6 ''^7 1 

17^*1 

31^ 

1 ;? 7 1 c 

4R M. 


1617 1 

240 2 

30 7 

1222-? 

A" "7 

-^7 

1 673 1 

36^7 

796 

^ 2 ? ? c 

72T 4 

-141 

1674 3 

4807 

?74 

7 2^1 q 

R4'-'4 

-18F 

16 3 8 3 

6 0 3 

746 

*» 2 7 c 

9A-V. 

-773 

1 A309 

720 6 

714 

^ ?qTC 

\r 

-767 

1 641 9 

8484 

1 RA 

f. 2S3<? 

I2'V’6 

-37^' 

' 64 5 c 

96 74 



1^7'7 


1 646 1 

10 808 


7 27C 1 

144^7 


^ 6479 

171’" 6 

96 

7 2 7 

1 '■? 


16^5 9 

1 3?0)7 

7r 

7 7R2 7 

1 3 ■■ 3 

-A r :j 

16 7 19 

144^^7 

48 

5 7^^ 7 

IR7 ’ ^ 

-67C 

1 681 c 

1 5 ^ ^ 8 

10 

1 2^7^ 

1971 

c 

1 68-^9 

16809 

-70 

7 2 ^<=> 7 

? 1 


1 68^1 

larr 9 

-66 

7 27^1 

7 1 0 1 

— R^C 

1703 1 

19719 

-10 7 

1317 7 

^73" 1 


1 7047 

20 45 9 

-14P 

1 ^*>77 

7 43| ? 

-^0^7 

1706 9 

21610 

-71 6 

1 3 2R7 

2^7^ 2 


1710 7 

77812 

-780 

7 30c q 

7SM 2 

— 1 /•* c r 

1711 5 

7A^>12 


1 343q 

27 A T 3 

-1 1 1 r 

1717 7 

25217 

-44 7 

1 3 /j c c 

7 8 3 1 4 

-1^61 

17207 

2641 3 

-51 3 

1 

3r ■’s 1 h 

-.1 2“9 

1 7367 

77 614 

-970 


31715 

-1 7 ^"4 

1 746 7 

28815 

-6^6 

7 3c 3 3 

3741 4 

-1 7Q7 

17527 

30013 

-7 37 

1 3^3*: 

3361"^ 

-1 33*^ 

17539 

3121 6 

-827 

1 3*31 « 

34 R1 7 

_1 R-7p 

17679 

32416 

-881 

1 37^!5 

3631 

-141 6 

17695 

33617 

-95 8 

1^875 

377-! 3 

-1 46 0 

17715 

3481 R 

-1030 

1 373 

3R 4’’ R 

-1 

1 7755 

36^'IB 

-10 96 

1 ^Q4 7 

396 70 

-1 546 

17763 

3721 8 

-1 166 

7. :s -^3 7 

408 7 1 

-1 

1 777 5 

3841 

-1234 

7 4 7 r 3 

47 \7 1 

- 1 6 2 7 

17885 

3962 3 

-1289 

1 7 3 

43 7/1 

- 1 A A 4 

1 8 1 5 

4^ 87 3 

-13?B 

*'-16 3 

4/ 4 2 2 

- 1 n 1 

18115 

42020 

-1 373 

1 1 7 1 

4567 3 

- 2 3 

1 8175 

43771 

-1621 

1 4 1 R 3 

46 8 7 3 

-1 ^76^ 

18 18 7 

444 7 2 

-1 476 

1 3 

48^20 

-*’ 7 ^'8 

1 8377 

45623 

-1 6-"'F 

1 447 3 

4922*^ 

PTTJ 

1 8 34 3 

4687 3 

-1562 



TABLE m.- TABULATED CONTROL POINT DATA 
FOR SAMPLE PROBLEM - Continued 


IX 

H . 

IZ 

IX 

IT 

IZ 








48-'^26 

-1 8 05 

2225 ® 

4562? 

-1 307 

1 2 

i^9225 

-1634 

22355 

46822 

-138? 

1 

506 26 

-1675 

22M 5 

8 7 2 3 

-1 45 6 


51626 

-1 714 

22429 

49 224 

526 

1 R ^ 

52876 

-1741 

22561 

5042 5 

>98 

1 B66 3 

540 27 

-17^7 

223:8 2 

51626 

-1663 

1R763 

55228 

-1779 

2260 3 

52826 

-172^ 

1RB23 

56479 

-1 8^*^ 

22 64 3 

54 27 

-1 786 

IRB3 B 

57629 

-1837 

22651 

55^27 

-1043 

^ 8075 

58829 

-1 340 

2266 3 

56-^27 

-1 8^7 

18091 

60P30 

-1 938 

777^: 

57627 

-1 086 

2001 1 

7 

216 

7 2 on 7 

58828 

-2022 

1 

1 20 0 

701 

2 3 *1 *> 3 

65 3 7 ) 

-7<-P2 

:>on«^q 

2430 

1 3 P 

A 3 

< 

“220 


363 7 

174 

2 ■'•.07 5 

1255 

^70t. 

2 018 1 

480 1 

156 

2421 5 

? V 3 f' 

— 004 

2031 1 

6 103 

1 4^ 

2/ 73 1 

3S^i 

— 1 07 

1 

72 13 

1 28 

7 /- 7 1 

/' 8 ' 7 

-182 

204-^1 

84^"^3 

ITT 

?62'31 

61 ^ 

«.i 6'- 

204Q3 

0604 

102 

74799 

72 ^'3 

-1 ^7 

20623 

108 38 

8P 

2^ 31 1 

8^v'>3 

-1*15 

20639 

120 )6 

69 

24391 

0 6 '\ f \ 

-7 A/. 

20689 

1320 6 

54 

2455 1 

1C 8 15 

-1 50 

208*=^ <; 

14 40 6 

38 

7 4 65 1 

J6 

-1 6 3 

20707 

156 07 

23 

2471 1 

137-^6 

— 1 

2r?i g 

^ 68'3q 

3 

2^*72 3 

144- 6 


2070 g 

180^' 9 

-20 

24 86 3 

1 

- 1 A 7 

209-^ q 

1920O 

-44 

2 4 P 7 c 

\(.oc ^ 

-161 

21089 

?04r*9 

-6 8 

7 4 80 q 

1 8 1 8 

-162 

21119 

21 61 0 

-07 

74939 

1070 8 

-1 71 

2 113 1 

2 2 811 

-12*'’' 

740/. 7 

7/1/. -7, 0 

-’PI 

?1 271 

24111 

-1 04 

74759 

2 1 6 1 1 

- 1 0 7 

21 28 7 

2521 1 

-1 PQ 

7 ^ 3 9 

2 281 T 

— 70 A 

21307 

26<'0. 2 

-221 

75109 

2 00 1 1 

-771 

213/^7 

2761 3 

-?6C 

2^2^ 9 

pr, 7^ 2 

- 2 35 

21385 

28814 

"3 ''"'6 

'’ 3 35 c 

264'’ 3 

-2 56 

21367 

300 14 

-30^‘ 

37 1 

27 51 3 

— 2^1 

2 1 46 7 

3121 5 

-*’96 

2 5 5 11 

28P ^ 3 

— pop 

21 6P 7 

32416 


2^527 

3 

— 33'’ 

21707 

33616 

-^.97 

25 56 1 

31 ?i 5 

-365 

21767 

3481 6 

-554 

2 5 5 81 

324 1 6 

-AOO 

21779 

36^17 

-6] P 

75909 

33617 

- AAP 

21 01 9 

3721 7 

-6 84 

7 5 65 7 

34817 


21938 

38419 

-786 

7 5 6 1 

36*117 

-^3A 

21 98 5 

39620 

-88^ 

25B49 

37218 

-506 

2190^ 

40830 

--976 

25047 

3B41 9 

—6'^ a 

22or 3 

4202"^ 

-1065 

2AOC 7 

3 96 1 9 

-723 

2201 5 

43221 

-1140 

2601 9 

4181 9 

—7 ^ P 

22095 

4442 2 

-1 230 

26159 

^.7121 

-850 



TABLE m.- TABULATED CONTROL POINT DATA 
FOR SAMPLE PROBLEM - Continued 



ir 

IZ 


lY 

IZ 


43231 

-935 

70021 

4081.9 

-765 

?A19*5 

44^ 22 

*-997 

3rr>8 7 

42^20 

-830 


45633 

-1-194 

301 87 

43221 

-91 r* 


468 2 3 

-11 87 

30247 

44421 

-989 


48'^?/. 

-1276 

39 25 5 

45621 

-1068 


49334 

-13^7 

393CC 

468 2 2 

-1153 

'^6 '**9 5 

5t’^474v 

-1477 

3041 5 

48)23 

-1236 

7 1 o 5 

^51674 

-1 A07 

39435 

49224 

-1326 


92337 

« 1 0 

3aa7 8 

5 7424 

-1414 

3 ^->66 7 

549 27 

-163^ 

7f 4P 3 

51625 

-1 rn2 

1 

55337 

-1 6 “^7 

3f 49 

52 82 6 

-1589 

?Jy^2 3 

56^^27 

7.36 

30 <57 C 

560 27 

-1 676 


57S28 

-1 ^93 

30 73 r 

55227 

-1 76A 

3 

58339 

~1 

30 83*^ 

56427 

-1 84 c 

7 A 09 1 

4, r. '>. •> 

-*» P96 

3 0 8 9*? 

576 2 8 

-1 930 

770C 3 

0 

-^27 

79007 

58829 

-2006 

7TQP 3 

12 n 

-Tr r. 

31 ''41 

6003 ) 

-2089 

2 31^ 3 

3^.00 

-5 8*^ 

7 ;?r A 3 

»**. 

-1467 

7 3 7*i o 

3 S " 1 

- A ^ £v 


12or 

-1381 

333^3 

4 3 ’ 7 

-S39 

33133 

24 or 

-1307 

233] 5 

61' >3 

— ^“OC? 

3313 1 

363 1 

-1 ?33 


73‘'"3 

-577 

3 314 7 

48-02 

-1177 


P-V'-- 7 

-5 3A 

^122 3 

64 6.2 

-1117 

2 8/' 9 1 

C6 

«/.0 9 

3338^ 

7202 

-1069 

3337 1 

19.817 


3 24 8 3 

84^3 

-I'^IS 

38*^39 

123^5 


32 54? 

9634 

-956 

70*91 

1330P 

^A\£ 

325*^5 

1083 5 

-R98 

^3 S3 1 

1.440 6 

- 30P 

3360P 

1230 5 

-859 

73^01 

15607 

-3P7' 

-271 1 

13205 

-8 OR 

7 0 3 0 1 

1 

-7TP 

32’'3 1 

1440 6 

-760 

7 P OJ” X 

1 P. 3 3 P 

-371 

7 7^71 

156^7 

-71 2 

339S3 

1 9?"o 

-363 

32^^ “^9 

16867 

-659 

7 Q 1 3 

2-‘4lo 

-?61 

? 2‘'91 

18 O'"'?- 

-6^c 

301 1 9 

21610 

-3^6 

33871 

192^^3 

-567 

7 9 1 3 9 

22316 

-7^3 

33^31 

20410 

-534 

3 91 9 

24113 

-351 

3313 1 

21 61 ^ 

-503 

301 8 7 

2531 1 

- 

33191 

22310 

-AP6 

391^9 

26413 

-358 

37 7 

24011 

-479 

30379 

27613 

-3“^i 

3334 3 

25212 

-47-7 

9 

28814 

-386 

73359 

26413 

-472 

39979 

3’>oi 

— 4'‘-r 

33379 

276^. 3 

—46 8 

7C^=-OQ 

312^ 5 

-6.1 A 

37 A1 0 

2881 3 

-466 

7 0 f 1 1 

3341 ^ 

-4 28 

3 3427 

30014 

-464 

79-7*. 1 

336^ 6 

-Arp 

3343 9 

31215 

-46^ 

707^,7 

3-i 3^ 6 

-AO 7 

3351 9 

326 16 

-4^6 

\) 

0 

CD 

36117 

-535 

33679 

336^ 6 

- 691 

39337 

3731 8 

— 89 

33779 

348] 6 

-508 

29335 

38419 

70 

338^ 9 

36017 

-54? 

3 9 3^* 7 

3961 9 

-7A3 

3385 T 

37218 

-582 



TABLE m.- TABULATED CONTROL POINT DATA 
FOR SAMPLE PROBLEM - Continued 



IL 

IZ 


lY 

IZ 

3^<591 

3H4] 8 

-624 

3T679 

36018 

- 53^0 

^5.007 

3961 8 

-698 

37759 

372^ 9 

-'^••12 

3^^^27 

4081 9 

-766 

3791 c 

3 841 9 

-677 

34^6 7 

42)21 

-832 

38.H9 

39619 

-732 

3407 5 

43221 

-897 

3 807 9 

40870 

-786 

3<^^087 

44421 

-98*=: 

38091 

42021 

-837 

34167 

45622 

-1067 

3P73 1 

43221 

-RO] 

34^27 

46823 

-1144 

38247 

444 ? 1 


3442 7 

48123 

-1277 

3R267 

4 362 3 

-1354 

34437 

49223 

-1 341 

3 8 397 

4682 3 

-1127 

34409 

5('42 3 

-lAAC 

38315 

48024 

-119“’ 

34639 

51624 

-1. 540 

38327 

40374 

-131"' 

346*^5 

52826 

-1 64r 

38 40 7 

5( 4 25 

-1 A1 8 

3467*^ 

5 A3 27 

-1739 

38567 

51636 

-.1 K1 c 

34^^1 

36727 

-1 RA7 

3 866 7 

528 3 6 

_i 632 

34*^2 3 

56^ 27 

-193r 

38727 

f 4327 

TKp 


57628 

-7040 

3 8 739 

6*^777 

-1 3 30 

3‘^^.RK 5 

S8829 

- 2 1 ^ 

38879 

56438 

- 23 1 r 

3497 5 

61029 

-7242 

38 89 5 

67629 

-?i 

3 6'^7 5 

n 

-2 5*=^. 4 

3801 r 

*^8830 

-3372 

36135 

120 0 

^ 7 A 7 r 

3 8 9 5 

.6i'0'’0 

-3301 

361.4. 7 

2A1^ 

—22 89 

Anor.R 

24 ’:1 


362P7 

36^^ 7 

-21 95 

40 31 5 

26 V 

-3' 1 3 

36^H3 

48 '*^3 

_2^T7 

4031 5 

4R‘. 2 

- 33^6 

36323 

63''3 

- 1 9 3 7 

AO 37 5 

63 '■'•2 

_ 7 n -? 7 

36 36 3 

72' ’ 3 

-1 Qro 

403 3 7 

73 0 3 

— -> P '5’7 

36371 

84'^‘4 

-1 68C 

4 f 7 

84' < 

- ? 7 3 p 

3638? 

OSi’i 5 

1 ^ 5 c. 

4054? 

96 8 

_ 7 c r. n 

3 6463 

1 -8^5 

A 7 P 

1563 

108" 5 

-3 3 Po 

3f 62 3 

1 2">.-5 

-1.414 

40 3 

1 3 0 '' 5 

— .3 ? 3 P 

36-^23 

1320 6 

-1 34^ 

40 61 1 

13306 

-3 "’73 

36-^8 3 

1 4407 

-1 271 

40A73 


-1 Q38 

36706* 

18S0 8 

-1 1 86 

607 0 

16A 7 

— 1 • 7 qo 

369*> 5 

16808 

-1 115 

AO 8 A 3 

148'^7 

693 

3695 1 

1 800 q 

-in 30 

4r 96 3 

1 83 *.R 

-.1 SPp 

36971 

1921 1 

-OK r 

41 02 3 

192« 9 

-1 AA6 

3731 1 

2041 0 

-882 

41 035 

2 041 0 

— 1 R R 0 

3731 Q 

2161 1 

- 82 3 

41 175 

21610 

-1.291 

3703 1 

228*' n 

-769 

4 119 1 

2 2BJ }. 

-1 211 

371 1 1 

2401 1 

-725 

41211 

24 V 2 

-1 137 

37 271 

2521 3 

-6 8 8 

41251 

2521 3 


37371 

26413 

-64? 

41259 

2641 3 

-OPT* 

37431 

2761 3 

-^97 

41 271 

27613 

-91 2 

37443 

28814 

— K P {^. 

413 51 

28814 

-847 

3758 ? 

3001 S 

-*^7C 

6 15? 1 

3C01 3 

-701 

37599 

3171 6 

-563 

41611 

31 2 : 6 

-“^44 

3761 9 

32416 


4167 1 

32416 

- 697 

37659 

33616 

-553 

416 8 3 

3361 6 

-68*^ 

37667 

34317 

-549 

41823 

34817 

-679 



TABLE m.- TABULATED CONTROL POINT DATA 
FOR SAMPLE PROBLEM - Concluded 




IZ 


IK 

IZ 

M 83 c 

3603 8 

— 70 c- 

461.5 9 

45623 

-3 118 


37218 

-7?? 

46239 

46824 

-lion 


38^1 B 

-“^66 

46399 

480 24 

-1263 

/^lor.7 

3963 9 

-ft^l 

46499 

49225 

-1 336 


A A 9^7 

-83r 

465*^ 9 

50 476 

-1402 

-^1 9<=>q 

4202 3 

-886 

4657 1 

51626 

-1 AA7 


43721 

— 0^6 

4671 1 

52827 

-164? 

A ^ 

444 2 2 

-1.0 2*^ 

46727 

540 27 

-1625 

K- 2 3 1 Q 

43623 


46747 

55?27 

T03 

1 

44 823 

-11 40 

46787 

56429 

-1778 

f 2471 

A 8 2 3 

-.1 2TA 

46^9*= 

5763 0 

-1°87 

A74P 7 

49273 

-1 288 

468 0 7 

^8830 

-2188 

42S7 7 

8f 42 4 

-1.766 

468 8 7 

6^06 30 

-237? 

42^4 7 

51626 

-1 437 

48 86 7 

16808 

-7549 

A 2 ^ 

52826 

-*5 53? 

49on 7 

IBO-^Q 

-329A 

42^67 

84-^27 

— 1 AA5 

A 002 3 

102’ 0 

-3019 

A2AA 7 

3^277 

-1 834 

4 904? 

204|0 

-277? 


54A2 8 

-1071 

4 2 0 8 3 

2161 0 

-2^-82 

A?c>77 

5 7629 

-21 21 

49^'' 0 1 

22811 

-2^07 

A?9A 7 

3 ft 3 2 Q 

-22 60 

491 0 3 

2401? 

-2214 

4?77<; 

6^r 29 

-7A1 

4018 3 

282'' 2 

-7082 

AAAI *: 

l”^ftr'3 

-3539 

49743 

26412 

-1 961 

2 7 


-331 3 

49/43 

27613 

-1 81 8 

44 76“’ 

1 3707 

-2132 

495^3 

28814 

-1 686 

AA.'Tp 3 

14/-..' 9 

-7927 

4951 5 

30015 

-1 55A 

A. 4 8 7 

156* 8 

-27 1 p; 

49f ^ 5 

31215 

-t A68 

44 0 4 3 

IAB'' 9 

- 2 c: 0 r 

496-^1 

3 2416 

366 

4A8 5 1 

1 ft o 

-27P P 

49601 

3361 7 

-1 274 

448 A 3 

1<7?1 ^ 

-2085 

4973 1 

34817 


4 4, ^4 3 


-1077 

4^*73 9 

36017 

-1152 

4*^.1 3 

7161 0 

-1 83^: 

4075 1 

3721 8 

-1133 

4P 7r 3 

2 28^.1 

-1736 

AC Q 3 1 

38418 

-■»!2.6 

4 3 7 4 3 

2471 3 

-1 6.30 

49OQ 1 

3962 j 

-1 ’.30 

C ? 7 8 

7 5 .7 ’ 3 

-1^37 

s r 0 0 1 

41801 

-11/7 

4=^41 8 

2641 3 

-.’» 470 

5015 1 

42021 

-1 |-^6 

^^^,43 1 

276 1 4 

33^. 

501 63 

43221 

-1 2^6 

A'^Af 1 

288^5 

-1 23T 

5 0 3r 3 

44422 

-1243 

4 3491 

3C'" 1 6 

— ’’176 

5 03’ g 

45623 

-12 67 

4 54 0 q 

3121 6 

-VO 89 

50339 

46823 

-12 84 

48^1 3 

3241 6 

-987 

50379 

48'‘23 

-1 7/^A 

1 

33617 

— 076 

50381 

40224 

-1 3/6 

4578 \ 

3481 A 

-9-4 

5 0 3 9 9 

50425 

-1 A3A 

45 85 1 

3601 9 

— 8 6'0 

5 0 479 

61626 

-1514 

A 8 0 3 1 

3721 9 

-Q77 

S n ,3 -a 9 

52826 

-1 ft 0/4. 

4^92 8 

3^419 

-01 1 

50739 

54^77 

-** 672 

AS 74 3 

3^62 6 


?=:07o c 

5522B 

-17AC 

4 5 778 

4-.;A21 

— 9 P4 

1 

8642 8 

-183ft 

4i^ 7c 9 

A?- 21 

-1012 

5 0 0 1 

57628 

-1025 

4 517 9 

43791 

-1 0A2 

5^'-967 

58828 

- 7 7 0 

46 3 A 7 

4462 2 

— 1 5-J"7 

98 7 

60 0 79 

-209] 



TABLE IV.- DEFINITIONS AND MATHEMATICAL EQUIVALENTS 
OF MAJOR FORTRAN VARIABLE NAMES^ 


Primary 

FORTRAN 

name 

Equivalent 

FORTRAN 

names 

Mathematical 

equivalent 

Definition 

A,B,C 


■^i’®i’^i 

Arrays containing coefficients 
of gradient vector 

ARG 


Hi,Z. 

Array containing real and 
predicted scalar variations 

BB 


Ni 

Array representing esti- 
mated gradient 

CLASS2 


Uj 

Array defining categories of 
discrepancies 

CONT 



Array of contour lines inter- 
secting side 1 

CONTA 



Array of contour lines inter- 
secting side 2 

CONTB 



Array of contovtr lines inter- 
secting side 3 

CONTC 



Array of contour lines inter- 
secting side 4 

CONTD 



Array of contour lines inter- 
secting side 5 

COUNT 



Array containing a count of 
discrepancies (frequency 
of occurrence) 

D 



Array for storing in memory 
all the control -point data 

DEPTH 



Array within which the 
2SPP + 1 by 2SPP + 1 
elements of the submatrix 
are stored 

E 



Array of vectorial compo- 
nents of a vector approxi- 
mating gradient 


GRAD 



Array containing components 




of a vector approximating 




the gradient 

I 

MAT, IE L 

Hi 

Depth-matrix array 


^Variables in table II are not listed except where mathematical equivalents are given. 



TABLE IV.- DEFINITIONS AND MATHEMATICAL EQUIVALENTS 
OF MAJOR FORTRAN VARIABLE NAMES^ - Continued 


Primary 

FORTRAN 

name 

Equivalent 

FORTRAN 

names 

Mathematical 

equivalent 

Definition 

lADDI 



Constant bias added to I 

ICODE 



Signifies end of data 

IDE PH 



Array for storing three rows 
of I matrix 

IDIST 



Distance of control data 
point from a common 
origin 

IFN 



Array containing location of 
input and output files for 
SORT2 

ISM 



An array containing informa- 
tion regarding input record 
structure and number of 
input files and of variables 
to be ordered 

DC 

X 


X- coordinate of control - 
point data 

lY 

Y 


Y- coordinate of control - 
point data 

IZ 

Z 


Z-coordinate of control- 
point data 

KEY 



Array containing location 
and t3q)e of variable to be 
ordered and direction of 
sort 

KORE 



Locations of the X- 

coordinates of each cor- 
ner of I matrix 

KORN 



Locations of the Y- 

coordinates of each cor- 
ner of I matrix 

KM 

JJ,MAT 


Number of rows in I matrix 

KMAX 


k’ 

Number of neighboring con- 
trol points 


^Variables in table II are not listed except where mathematical equivalents are given. 
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TABLE IV.- DEFINITIONS AND MATHEMATICAL EQUIVALENTS 
OF MAJOR FORTRAN VARIABLE NAMES^ - Continued 


Primary 

FORTRAN 

name 


Equivalent 

FORTRAN 

names 


Mathematical 

equivalent 


Definition 


KMAXl k 

LM KK,MIT 


Number of neighboring con- 
trol points 

Number of columns in 
I matrix 


M 

MAXWRD2 

MAXWRD3 

MINWRD2 


Number of unfilled I elements 
Maximum lY 
Maximum DC 
Minimum lY 


MINWRD3 Minimum DC 


Locations of X-intercepts 



with side 1 


POSIY 

Locations of Y 
with side 1 

intercepts 

POS2X 

Locations of X 
with side 2 

intercepts 

POS2Y 

Locations of Y 
with side 2 

intercepts 

POS3X 

Locations of X 
with side 3 

intercepts 

POS3Y 

Locations of Y 
with side 3 

intercepts 

POS4X 

Locations of X 
with side 4 

intercepts 

POS4Y 

Locations of Y 
with side 4 

intercepts 

POS5X 

Locations of X 
with side 5 

intercepts 

POS5Y 

Locations of Y 
with side 5 

intercepts 


Radial distance between 
^ two points 

^Variables in table n are not listed except where mathematical equivalents are given. 
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TABLE IV.- DEFINITIONS AND MATHEMATICAL EQUIVALENTS 
OF MAJOR FORTRAN VARIABLE NAMES^ - Concluded 


Primary 

FORTRAN 

name 

Equivalent 

FORTRAN 

names 

Mathematical 

equivalents 

Definition 

RLAST 



Inner radius of ring 

RMAX 



Maximum distance of data 
point from arbitrary 
origin 

RSAVE 



Outer radius of ring 

SP 



Distance between elements 
of I at plotting scale 

SQX 



X- coordinate of submatrix 
center 

SQY 



Y- coordinate of submatrix 
center 

W 



Weight 

XMAX 



X-dimension of I 

YMAX 



Y-dimension of I 


^Variables in table II are not listed except where mathematical equivalents are given. 
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